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Abstract 


The  UTIAS  10  cm  x 18  cm  Hypervelocity  Shock'  -Tube  has  ^ been  used  In 
recent  years  to  study  ionizing  shock  structures,  flat-plate  and  slde-vall 
boundary  layers  and  nonstationary  oblique  shock-wave  diffractions  over 
compression  comers.  These  phenomena  were  recorded  using  a 23-cm  dia  field 
of  view  Mach-Zehnder  interferometer  equipped  with  a giant-pulse  dual-frequency 
ruby  laser. 

In  order  to  extract  Ithe  maximum  amo\mt  of  data  in  these  complex  flows 
a digital  evaluation  method  was  employed.  For  this  technique  a new  approach 
to  the  theory  of  Interference  was  developed.  In  this  approach  the  spatial 
coordinates  (x,y)  of  the  various  lines  of  Interference  (fringes)  on  the  inter- 
ferograms  are  put  into  digital  form,  thereby  making  a computer  analysis 
possible . 

The  experimental  technique  and  insbrimientation  associated  with  the 
various  measurements  Involved  in  research)- r using  the  UTIAS  10  cm  x l8  cm 
Hypervelocity  Shock  Tube  are  described.  Finally,  the  maximum  possible 
absolute  and  relative  errors  associated  with  these  measurements  are  cidculated 
for  four  different  gases  (argon,  krypton,  oxygen  and  nitrogen)  which  are 
usually  used  as  test  gases. 
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Notation 

a entrance  position  to  test  section 

a speed  of  sound 

Ar  argsn  atom 

argon  Ion 

b exit  position  from  test  section 

c^  speed  of  light  in  vacuum 

d unsupported  diameter  of  the  diaphragm  » 18.1  cm 

e electron  charge 

E(0)  absolute  error  In  measuring  quantity  (ji 

H the  height  difference  of  the  oil  In  the  manometer 

2 

g standard  acceleration  of  gravlty(*  980.66^  cn/aec  ) 

h resldvial  thickness  of  the  dit^hragm  after  scribing  (Xq.  1.2) 

k,  arbitrary  constant 

K constant  (Eq.  1.2) 

Gladstone>Dale  constant  of  species  J 

Kr  krypton 

!•  ' 

.1  geometrical  ItoXthctravellsd'oby  light  .‘section 

L geometrical  distance  across  test  section 

m^  mass  of  species  J 

M continuous  fringe  nuniser 

M'  integer  fringe  number 

M incident  shocks  wave  Mach  nuaiser 

s 

n index  of  refraction 

n constant  (Eq.  1.2) 

n^  nuniser  density  of  species  j 

! 

N continuous  variable 


i 
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Notation  - continued 

N • integer  variable 

N2  nitrogen  molecule 

^2  oxygen  molecule 

P pressure 

P^urat  pressure  of  the  diaphragm 

P^  initial  pressure  ahead  of  the  diaphragm 

^hurst 

R gas  constant 

S fringe  shift 

5 distance  hetveen  two  shock  detectors 

t time 

t total  diaphragm  thickness  (£%.  3i.2} 

T temperature  in  °K 

initial  temperature 

laboratory  temperature 

velocity  of  incident  shock  -wave 
a 

X spatial  coordinate 

X collision  partner 

y spatial  coordinate 

z spatial  coordinate 

0!  degree  of  dissociation 

^ induced  electric  dipole  polarizability  of  the  gas 

7 specific  heats  ratio  "(Cp/Cy) 

6 first  variation  of  an  integral 

A0  change  in  property  0 

e angle  of  rotation  of  second  beam  splitter  in  interferometer 
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Notation  - continued 


9^  wedge  angle 
X light  source  wavelength 

X^  light  source  wavelength  in  vacuum 

P density 

P^il  density  of  the  oil  in  the  manometer 

T optical  path  length 

1'*  optical  path  length  inside  the  test  section 

light  source  frequency 

plasma  frequency 

X degree  of  ionization 

Subscripts 
a atom 

e electron 

E ionization  equilibrium  conditions 

1 ion 

m molecule 

p test  section  arm  of  the  interferometer 

q compensation  chamber  arm  of  the  interferometer 

r reference  point 

o "no-flow”  conditions 
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1.  INTRODUCTION 


The  UTIi^  10  cm  X l8  cm  Hypervelocity  Shook  Tube  was  built  in  19^5 
and  reported  in  detail  (design,  instrumentation  and  performance)  by  Boyer 
(Ref.  1}  upon  completion  of  the  construction  and  calibration.  In  addition 
Bristow  (Ref.  2)  reported  on  further  modifications  and  iaqprovements . 

Since  then,  further  modifications  and  changes  were  made  by  various 
researchers  who  have  used  the  facility  for  investigating  dissociating 
and  ionizing  gases  and  various  shock-wave  phenomena. 

Since  the  original  design  was  in  British  units,  w^ilch  were  replaced 
recently  by  the  metric  system,  and  since  the  various  modifications  are 
described  in  various  UTIAS  publications,  it  was  felt  necessary  to  combine 
them  Sbll  together  in  an  updated  document  that  will  describe  the  present  status 
or  the  shock-tube  as  well  as  its  maintenance  and  operation  procedures. 

The  investigated  phenomena  are  usually  recorded  using  a 23  cm  dla 
field  of  view  Mach-Zehnder  interferometer  equipped  with  a dual-frequency 
glaht-pulse  laser.  Consequently,  the  basics  of  two-wavelengthi'  inter- 
ferometry are  reviewed,  and  a new  and  very  accurate  interferometric  data- 
evaluatlon  method  is  presented. 

2.  UTIAS  10  cm  x l6  cm  HYPERVELOCITY  SHOCK  TUBS 

The  UTIAS  combustion  driven  hypersonic  shock  tube  is  a large , well- 
constructed  facility t approximately  IS  meters  in  length,  designed  to 
simulate  high  enthalpy,  real  gas  flows  typical  of  hypersonic  flight.  !Ihe 
coBibustion  driver  is  capable  of  generating  shook  speeds  in  excess  of  7 
kilometers  per  second.  A detailed  account  of  the  component  design  and 
construction  features  has  been  given  by  Boyer  (Ref.  l),  while  improvements 
and  modifications  to  the  system  have  since  been  reported  by  Bristow  (Ref.  2) . 

A schematic  overview  of  the  shock  tube  facility  in  its  present  form 
is  shown  in  Figure  1.  In  describing  the  major  components  and  operation,  an 
attempt  is  made,  where  possible,  to  convert  the  dimensions  of  feet  and  Inches 
from  the  original  design  to  more  conventional  metric  units. 

2.1  Shock  Tube  Driver 


The  chrome-molybdenum  steel  cylindrical  driver  has  an  internal 
diameter  of  13.24  cm  (6”)  wlth^a  wall  thickness  of  10.16  cm  (4”),  and  can 
withstand  pressures  up  to  11,000  psl  with  a safety  factor  of  6.  A driver 
length  of  1,422  m (36")  was  considered  adequate  and  economical  for  the 
experiments,  although  other  lengths  are  available.  Driver  pressures  up  to 
about  1000  psi  are  usually  acquired  directly  from  gas  bottles  while  higher 
pressures  are  obtained  by  the  combustion  of  a stoichiometric  mixture  of 
hydrogen  and  oxygen,  diluted  with  helium.  The  optimum  mixture  used  for 
smooth  combustion  and  r<*lat  vely  fast  rise  time  was  determined  by  Bristow 
(Ref.  2),  which  consisted  of  oxyc^n,  20%  hydrogen  (includes  3%  used 
as  a partial  diluent)  and  ’12 A%  helium  by  pressure.  As  indicated  by  Benoit 
(Ref.  3),  the  pressure  should  increase  by  a factor  of  about  8 upon  combustion, 
although  experiments  have  shown  this  factor  to  be  around  6.7.  Consequently, 
a conservative  estimate  of  6.1  was  used,  giving  an  insurance  factor  of  lOit 
when  determining  the  correct  mixture  pressure  for  a particular  diaphragm 
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bursting  pressure. 

The  combustion  is  achieved  by  heating  a long,  0.38  mm  (0.015”) 
diameter  tungsten  vire,  suspended  over  the  length  of  the  driver  along  the 
centre-line.  To  obtain  the  necessary  wire  teagjerature , a 4.5  Hf  capacitor, 
maintained  at  13.4  Kv,  is  discharged  through  a triggered  gas-type  spark  gap 
and  the  wire  to  ground.  Poor  Ignition  is  minimized  by  adequately  Insulating 
the  high  voltage  connection  through^Ahe  driver  wall  with  a continuous  4 mm 
thick,  teflon  sheath,  and  ensuring  a bright  yellow  glow  of  the  wire  under 
atmospheric  testing. 

To  obtain  excellent  repeatability  and  control  of  the  combustion 
pressure,  the  driver  is  evacuated  to  approximately  1 torr  prior  to  each 
experiment,  using  a Cenco,  Byvac  l4,  mechanical  vacuum  pump.  %e  correct 
amount  of  oxygen  for  the  final  mixture  is  then  admitted  and  critically 
monitored  to  + 0.1  psi  using  a Heise,  H17507C,  dial  pressure  gauge.  Helium 
is  then  added*to  reach  46j!  of  the  mixture  pressure,  hydrogen  to  reach  66}{, 
and  finally  helium  again  to  complete  the  mixture.  By  using  this  sequence 
and  admitting  all  gases  through  a mixing  tube  with  small  orifices  every  15  cm 
along  the  complete  length  of  the  driver,  excellent  combustion  performance 
la  ensured. 

Diaphragms 

The  diaphragm  clamp  between  the  driver  and  driven  sections  is 
designed  for  diaphragm  discs,  35.4  cm  in  diameter  (13.875”) » sad  leaves 
an  unsupported  diameter  of  l8:.l  cm  (7.125”).  0-rings  seal  against  the 
supported  annulus  of  the  diaphragm  on  each  side.  Two  types  of  diaphragm 
can  be  used*  for  non-combustion  driver  pressures  less  than  1000  psi 
(M  < 10)  mylar  discs  are  out  from  sheets  ranging  in  thickness  from  0.001” 
to  0.050”.  For  the  unsupported  diameter  of.  l6.1  cm,  mylar  was  found 
experimentally  to  have  a bursting  pressure,  varying  linearly  with  thickness, 
as 


^urst  " i 


(1.1) 


\niere  t " mylar  thickness  in  Inches.  Consequently,  desired  bursting 
pressures  can  be  achieved  merely  by  combining  discs  of  various  thickness . 

It  should  be  noted  that  the  sylar  diaphragms  burst  without  shattering, 
leaving  the  tube  free  from  debris.  The  longer  formation  distance  caused 
by  a slower  opening  is  not  a problem  in  this  long  facility. 

The  stronger  shock  waves  are  generated  using  annealed,  stainless 
steel,  type  304  (SS  304)  diaphragms  with  a 2-^  finish  to  provide  adequate^ 
vacuum  sealing.  These  diaphragms  are  machine-scribed  in  the  form  of  a 90° 
cross  on  the  driven  side  to  enable  petaling  rather  than  shattering  eiien 
burst  by  the  combustion  pressure.  General,  thicknesses  in  the  nangsi  0;062”  to 
0.182”  are  sufficient  to  cover  the  range  of  necessary  bursting  pressures 
for  10  < M < l8  by  machining  the  scribe  depth  appropriately.  For  each 

mm  9 

diaphragm,  the  bursting  pressure  can  be  predicted  from  the  calibrations 
of  Bristow  (Ref,  2)  as; 


Wst  "'  ( ^ ) (I)  »“ 

where 

t - total  diaphragm  thickness 

h - residual  thickness  of  the  diaphragm  after  scribing, 

d » unsupported  diameter  of  the  dla:^ragm  ■ l6.1  cm, 

^ - ultimate  tensile  strength  of  the  annealed  material, 
which  for  SS  30l»  is  65,000  pal. 

The  remaining  parameters  were  found  experimentally  from  the  calibration  as: 
n ■ 2.2  + 0.1 

K = 1.0l"for  t a 0.172"  and  d * l8.1  cm 

K ■ 1.10  for  t - 0.109"  and  d ■ l8.1  cm 

■ 1.06  for  t ■ 0.062"  and  d ■ l6.1  cm 

To  ensure  proper  petaling,  the  ratio  of  residual  thickness  to  totislvthlokness 

Is  kept  between  0.55  j<h/t  <0.75.  The  warlatlon  In  total  thickness  over 
each  diaphragm  is  kept  belcw  Q.004"  while  the  residual  thickness  is  machined 
to  an  accuracy  of  0.0005"  using  a vacuum  chuck.  These  rigid  tolerances  are 
maintained  both  for  safety  and  for  repeatability  of  experiments . Machining 
details  have  been  outlined  by  Boyer  (Bef.  l). 

For  optimum  performance,  the  diaphragm  is  chosen  to  have  a bursting 
pressure  equal  to  the  driver  pressure  P]^  required  to  drive  a shook  wave  at 
specified  Mach  number  N into  the  test  gas  at  specified  initial  pressure  p. . 
Hydrogen,  helium  and  ca»on  dioxide  are  usually  used  as  the  driver  gases 
for  obtaining  shock  waves  In  the  range  2 < <10  with  mylar  diaphragms . 

In  this  ease,  the  bursting  pressure  Is  determined  from  standard  curves 
relating  the  required  pressure  ratio  P]^/P2^>  to  the  desired  Mach  number 

(Ref.  1:)  with  a correction  factor  of  about  2 used  to  account  for  non-ideal 
effects.  A recommendation  of  driver/driven  gas  combinations  and  the  required 
pressure  ratio  pj^^  for  Incident  shook  waves  in  the  range  2 ^ ^ 8 for 

oxygen,  nitrogen  and  argon  aa  test  gases  is  given  In  Table  1.  For  combustion 
driver  experiments  into  argon  Bristow  (Ref.  2)  developed  the  following 
emplriosl  e:qpresslon 


in  ■ 0.422  Mg  + U.02  10  < M^  < 2k 

A different  expzwsslon  was  claimed  by  Whitten: 


1350  (M  -10)^,  for  12  < M < 18 
s — s — 


Unfortunately,  a slnllar  expression  for  krypbon  was  not  developed,  due  to  the 
lack  of  sufficient  nutaiber  of  experiments.  Hovever,  Tang  (Ref.  5)  established 
a empirical  relation  for  M « 15 


(0.422  + 2.92) 

s 


above  mentioned  expression  is  P.  ^ as  calculated  from 

burst 


Shock-Tube  Driven  Section 

As  shown  in  Figure  1 the  driven  section,  made  from  SPS-245  tool 
steel,  has  a total  length  of.  15.24  m (50'),  including  a 0.457  n (1.5 ') 
transition  section  from  circular  to  rectan^ar  cross-section,  and  a 
specially  designed  1.2192  m (4')  test  section  located  13.56  m (44.5*)  from 
the  dia^ragm  stetlon.  It  is  tarminatad  by  a large  dump  tank  in  the  form 
of  a vertical  cylinder  with  hemispherical  caps  whose  internal  diameter  is 
0.9144  m (3'),  overall  height  is  *.93  m (6.5')  and  whosa  volume  is  1 m^ 

(35.4  ft3).  The  rectangular  channel  cross-section as  indicated,  has  a 
horizontal  width  of  10.16  cm  (4")  and  a vertical  height  of  17 .73  cm  (7**). 

The  internal  surfaces  were  originally  ground  to  a 0;.8  micron  (32  microinch) 
finish  and  plated  with  a 51  micron  ( .002”)  layer  of  chromium.  Instrumentation 
ports,  lettered  A,B,  C,,.. (Figure  1)  are  located  at  1.2192  meter  (4*) 
intervals  along  the  drlvea  section  allowing  for  installation  of  necessary 
gauges. 


Note  that  Pi,  in  the 


The  test  section  contains  6 instrumentation  ports,  and  incorporatea 
two  interferometric  quality  circular  windows,  basically  20.32  cm  (8”) 
in  diameter  end  9.53  cm  (3.75")  thick.  These  windows  vere  made  from 
borosilicate  crown  (BK-7)  optical  quality  glass,  with  surfaces  ground 
and  polished  flat  to  within  a wedge  angle  less  than  2 seconds  of  arc,  and 
specified  to  have  a maximum  peak  to  peak  variation  in  transmitted  wavefront 
over  t^  entire  diameter  of  less  than  one-quarter  of  a wavelength  at 
6328  A°.  These  rather  large  windows  were  very  carefully  mounted  in  frames, 
with  flat  neophrene  rubber  vacuum  gaskets  and  25  micron  removable  shims  used 
to  prevent  any  metal  to  glass  contact.  The  sensitive  mounting  is  iaqportant 
with  this  size  of  window  in  preventing  az:y  wdue  stress  concentration  and 
cracking  from  either  the  shock  loading  or  the  initial  mounting  itself,  and 
essential  from  the  standpoint  of  safety,  economics  and  interferometric 
quality. 


Surface  burns  and  deposits  that  mi^t  develop  on  the  windows, 
particularly  near  any  perpendicular  metal  surfaces  (test  models)  are  often 
removed  using  a dilute  solution  of  hydrochloric  acid,  and  were,  if  necessary, 
moved  from  a particular  region  of  Interest  by  a l8o°  rotation  of  the 
rectangular  window  frames. 
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The  entire  driven  section  and  dun^  tank,  coomrising  a volume  of 
1.28  m-’  (45.13  ft^),  are  evacuated  prior  to  each  experiment  to  maintain 
the  highest  possible  purity  of  the  test  gas  under  consideration.  The  pum]p~ 
Ing  system  Includes,  in  series,  a Kinney  type  KS-«47  mechanical  pump,  a 
Heraeus  Roots,  type  RG-350,  blover  (forepump)  and  a CVC,  type  PMCS-6B, 
six-inch  diameter  diffusion  pump.  The  ultimate  vacuum  after  a pumpdown 
period  of  about  three  days  was  3 x 10“^  torr,  with  a leak-outgasslng  rate 
of  2 X 10”5  torr  per  minute.  The  slightly  longer  pumping  time  than  that 
reported  by  Bristow  (Ref.  2)  can  be  attributed  to  a modification  of  the 
coupling  between  the  vacuum  punp  system  and  the  shock  tube.  Previously,  the 
vacuum  manifold  had  been  c(mnected  to  the  tube  by  three  pipes,  5.08  cm  (2") 
in  diameter,  to  ports  situated  between  1.5  m and  2.7  m from  the  diaphragm 
station.  A high  pressiire  and  vacuum  service  ball  valve  was  located  on  each 
pipe,  approximately  21.3  cm  from  the  Interior  wall  of  the  shock  tube, 
allowing  a total  cavity  volume  of  1240  cm^  to  the  tube  at  all  times.  In  view 
of  the  instabilities  observed  by  several  experimenters  in  this  facility  for 
strong  ionizing  shock  waves  in  argon,  it  was  useful  to  eliminate  any  possible 
contributing  factors,  no  matter  how  remote.  Conseqtuently,  only  one  port 
was  retained,  with  the  ball  valve  modified  to  lie  about  5 cm  from  the 
interior  wall.  Ihe  remaining  ports  were  blanked  off  leaving  an  exposed 
cavity  volume  of  103  em^,  or  one-tenth  of  the  original  cavity.  Although 
no  noticeable  differences  were  observed  with  the  shock  stability  problem, 
the  single  connection  was  retained  as  the  more  appropriate  system,  particu- 
larly since  the  conductance  loss  did  not  seriously  affect  the  pumpdown  time. 


Test  Q^s  Conditions 

To  ensure  a reasonable  degree  of  cleanliness  before  each  experiment, 
a general  procedure  is  followed  in  wliich  the  driven  section  is  purged  with 
high  pressure  air,  the  test  section  windows  (in  their  frames)  are  removed 
and  cleaned  with  optical  solution,  and  the  test  section  walls  are  scrubbed 
in  succession  with  trichloroethylene,  acetone  and  finally  methyl  alcohol. 

(On  occasion,  the  entire  driven  section  is  cleaned  in  this  manner).  As 
mentioned,  the  driven  section  is  then  evacuated  over  a period  of  usually 
three  day**  Before  the  Introduction  of  the  test  gas,  the  pressure  is  monitored 
by  an  ionization  gauge  (Associated  Electrical  Industries,  Model  VC-10) 
located  at  an  instrumentation  port,  2.4  m from  the  diaphragm  station.  The 
pressure  in  the  vacuum  manifold  can  also  be  checked  with  a CVC  type  OPH-IOOA 
discharge  gauge.  A pressure  of  less  than  5 x 10’^  torr  is  generally  required 
before  proceeding  with  an  experiment.  A mass  speetrometric  investigation 
of  the  driven  section  under  these  conditions  indicated  the  pressure  to  be 
made  up  mainly  of  water  vapour  (70K),  nitrogen  (20){)  and  oxygen  (5^)*  In 
addition,  the  leak  rate  after  feoLation  of  the  vacuum  pumps  was  found  to 
be  approximately  2 x 10"5  torr /minute,  with  the  increase  in  pressure  due 
almost  entirely  to  air  (Eg  and  Og). 

In  the  experimental  procedure,  the  test  gas  is  admitted  Immediately 
after  Isolation  of  the  vacuum  puaps  to  minimise  effects  of  outgasslng,  and 
the  shock  tube  is  "fired"  approximately  2 to  5 minutes  following  this. 

Hence,  a maximum  impurity  partial  pressure  of  about  2 x lO*^  torr  is  ensured, 
\dilch,  for  a relatively  low  initial  pressure  of  5 torr,  amounts  to  an 
impurity  level  of  40  ppm.  As  the  test  gas  is  admitted  to  the  driven  section, 
the  pressure  is  monitored  at  the  vacuum  manifold  with  a Wallace  and  Tlernan, 
type  FA  l60  (0-50  torr)  dial  gauge.  A more  accurate  measurement  of  the 
initial  pre-shock  pressure  is  then  made  using  an  oil  manometer  (0.2-40  torr). 
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described  by  Bristow  (Ref.  2).  Ibis  device  consists  of  a $0  ea  Icmg*  1 cm 
inside<4iaaeter,  glass  U>tube  with,  at  one  end,  a bulb  reservoir  containii^ 
Dow  Coming,  DC  704  diffusion  pungp  oil,  and  at  the  other,  a flexible  connec- 
tion to  the  system.  The  manometer  is  vacuum  pumped  before  use,  keeping  the 
oil  in  the  reservoir  by  a 120°  rotation  of  the  tube.  In  the  i;^rlght  position, 
pressure  can  then  be  determined  simply  by  a sensitive  measurement  of  the 
height  difference  H,  between  oil  levels  in  each  arm,  i.e.,  p « P ..  gH  ^es/ 
cm^,  where  g » 980.655  om/aec^  (standard  acceleration  of  gravity)”^ 

pQii  - 1.069  + 9.5  X 10-^  (25°C-T)g/em3,  is  the  density  of  DC  704 
oil,  adjusted  about  c.5°C  for  theimal  expansion. 

T is  the  oil  temperature  in  the  manometer. 

In  units  of  torr,  p ■ P^^^H/13«59?1  where  H is  meuured  in  mm. 

The  temperature  of  the  driven  section  is  measured  to  an  accuracy  of 
+ 0>.1°C  prior  to  each  experiment  using  a standard  mercury  bulb  thermometer, 
rnserted  in  a 5 cm  deep,  6 mm  diameter  oil-filled  port  in  the  upper  wall 
of  the  ahobk  tube  walls,  this  measure^n't  is  taken  as  the  initial  pre-shock 
temperature  of  the  test  gas. 

2.2  Basic  Instrumentation 


Fast-response  electronic  equipment  is  used  to  monitor  the  shock-tube 
performance  in  each  experiment.  Ihls  includes  a shock  detection  system  for 
both  accurate  measurement  of  shock  velocity  and  microsecond  control  of  the 
diagnostic  equipment,  and  the  recording  -of  pressure  history  in  boLb  the 
driver  and  driven  sections. 

Shock  Detection 


In  each  experiment,  the  arrival  of  the  shock  front  at  several 
locations  along  the  driven  section  is  detected  using  Atlantic  Research, 
type  LD-25,  piezoelectric  pressure  transducers,  flush  mounted  in  selected 
instrumentation  ports.  These  gauges  have  a sensitivity  of  approximately 
0.15  V/psi,  a rise  time  less  than  1 microsecond,  and,  as  shock  detectors, 
require  a minimum  amount  of  associated  electronic  equipment . The  6 mm 
diameter  surface  of  each  gauge  is  given  additional  protection  with  a thin 
coating  of  silicone  rubber  (G£  RTV). 

The  transducer  output  is  connected  directly  to  a power  preamplifier*, 
situated  imediately  adjacent  to  the  gauge,  which  consists  simply  of  a 
field-effect  transistor  for  high-impedance  input  as  a source  follower, 
and  a me-stage  emitter-follower.  Although  the  voltage  gain  is  less  than 
unity,  the  power  gain  and  low  Impedance  output  enables  the  subsequent 
transi^ssion  through  long  lengths  of  cable  with  no  oppreciable  signal  loss. 

If  necessary,  the  signal  is  increased  in  voltage  at  v.he  main  electronics 
console  by  a small , two-stage , resistance-coupled  amplifier*  which  uses  a 
'*diode-catching*'  technique  to  give  a 20-volt  output  for  all  amplified  signals 
greater  than  a pre-set  level,  and  0 volts  for  any  lover  signals.  Correct 

*'  dourtesy  of  Mr.  A.iPerrin,  UTIAS 
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selection  of  this  level  before  each  experiment  effectively  discriminates 
against  anticipated  input  noise,  and  the  20-volt  output  signal  is  suitable, 
for  example,  in  triggering  conventional  electronic  counters.  The  overall 
response  tiiM  of  the  system  is  less  than  1 microsecond.  The  transducer, 
along  with  relatively  minor  electronic  equipment,  provides  an  economical, 
reliable  and  accurate  method  of  shock  detection  with  proven  durability.  The 
system  schematic  is  included  for  general  information  in  Figure  2, 

The  detectors  were  used  both  in  the  measurement  of  shock  velocity 
and  in  the  accurate  timing  of  the  laser  light  sourpe  for  interferogram 
exposure  discussed  in  Section  2.4. 

Shock  Velocity  Measttrements 

The  incident  shock  velocity  is  monitored  using  an  x-t  arrangement 
(a  typical  arrangement  is  dsowi  in  Fig.  3).  Three  Hewlett-Packard  (Type  3734A) 
and  two  Racal  (Type  SA.4^)  digital  counters  are  connected  to  a common,  external, 
one  megacycle  oscillator  giving  all  counters  a uniform  time  base  with  a one 
microsecond  resolution.  When  the  shock- wave  arrives  at  station  D,  (6.^  iii 
downstream  of  the  diaphragm  and  5.?  m upstream  of  the  test  section)  time 
counters  1,2,3 ,4  and  5 are  triggered  simultaneously.  These  five  counters 
are  then  stopped  in  succession  as  the  shock-wave  passes  subsequent  detectors 
at  stations  F,Q,H,I  and  J.  An  additional  Hewlett-Packard  (Type  ?325  A) 
counter  (No. 7,  Fig. 3)  is  used  between  stations  I and  J to  indicate  the  time 
taken  for  the  shock-wave  to  travel  past  the  test  section  and  the  test  model* 
Since  the  shock-wave  velocity  changes  somewhat  due  to  diffraction  over  test 
models  (Refs.  6 and  7),  the  measurement  of  counter  No .7  is  used  only  as  a 
monitor  of  the  shock  location. 

The  control  of  the  laser  light  source  operation  is  also  Included  in 
Fig.  3.  The  laser-flashlamp  capacitor  bank  is  usually  triggered  by  the 
shock  arrival  at  station  F,  suitably  delayed  in  a Tetronlx  type  55?  oscilloscope 
to  get  900  microseconds  of  energy-pumping-time  by  the  time  the  Pockels-cell 
shutter  opens..  The  Pockels  cell  Q-swltch  is  always  triggered  from  station 
I,  delayed  correctly  In  a pulser  delay  unit  to  take  the  Interferogram  at  a 
desired  time  (or  location  of  the  shock  wave  wrt  station  I), 

To  monitor  the  actual  flashlamp  pumping  time  in  each  experiment, 
synchronization  pulses  from  both  the  flashlamp  capacitor  bank  and  the 
Pockels-cell  pulser  unit  are  used  to  'itart  and  stop,  respectively,  a UTIAS 
microsecond  counter  (Fig.  3,  Counter  No. 6). 

As  an  additional  check  on  shock-tube  perfomance  in  each  experiment, 
a Kistler  type  601-B,  piezoelectric  pressure  transducer  is  used  to  monitor 
the  pressure  variation  with  time  behind  the  incident  shock  wave.  The 
transducer  output  is  taken  directly  to  a charge  amplifier  (Kistler,  Model 
504}  with  a final  signal  displacement  on  the  designated  oscilloscope  (Fig.  3). 

With  the  above  configuration,  an  x-t  plot  for  a typical  experiment 
is  shown  in  Fig.  4,  with  the  (0,0)  reference  taken  as  station  D. 

The  shock  velocity  along  the  driven  section  is  determined  in  two 
ways . The  average  velocity  between  stations  is  immediately  obtained  by 
finding  the  slope  of  straight  line  segments  Joining  successive  points  (i.e., 
a simple  division  of  the  distance  between  stations  by  the  difference  in 
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the  recorded  shock  aurrival  times  at  respective  locations) . The  error  in 
this  measurement  depends  on  the  maximum  spatial  uncertainty  of  + 6 mm 
arising  from  the  finite  size  of  the  shock  detectors,  and  the  + 1 ^pec 
uncertainty  in  time  neeordlhgj;.  For  the  do  jest  inter-station  distance  of 
0.6l  m,  and  a time  measurement  of  110  (isec  for  a typical  shock  speed  of 
5500  m/sec,  the  error  in  measured  average  shock  velocity  would  be,  at 
most,  2%  (see  subsequent  discussion  in  Section  2.3). 

A smooth  curve  is  also  fitted  to  the  set  of  (x,t)  points  by  a method 
of  least-squares  using  Chebychev  polynomials.  The  best  fit  of  lowest  order 
satisfying  the  uncertainties  in  each  point  is  determined  numerically,  and 
the  result , in  most  cases , is  a parabolic  or  second  order  fit . The  instan- 
taneous shock-wave  velocity  coUld  then  be  evaluated  from  the  derivative 
at  anytime  or  position.  Unfortunately,  the  differentiation  of  such  a fit 
always  involves  a loss  of  accuracy  which  is  somewhat  difficult  to  quantify, 
particularly  with  only  five  data  points  involved.  However,  with  the  smooth 
fit  satisfying  the  uncertainties  in  all  points,  it  is  felt  that  this  method 
is  slightly  better  than  the  calciL}.ation  of  average  velocities,  with  the 
accuracy  estimated  to  be  under  1%  within  the  bounds  of  the  data. 

Past  examination  of  the  Instantaneous  and  average  velocities  along 
the  driven  section  indicated  that  the  shock  speed  was  decreasing  or 
attenuating  as  expected.  For  example,  over  the  ler^h  of  5.49  m (between 
stations  D and  l),  the  maximum-  observed  change  in  wave  speed  was  found  to 
be  5/f  (often,  it  is  much  less).  In  other  words,  the  maximum  attenuation  is 
less  than  1)(  per  meter,  and  Just  above  the  uncertainties  in  the  velocity 
measurements  themselves.  As  a r«.  ilt,  for  each  experiment,  the  shock  velocity 
used  in  ensuing  cAlculations  is  ta^  'u  ' i the  mean  value  of  the  instantaneous 
velocity  determined  at  station  I imnuaiately  ahead  of  the  test  section, 
and  the  average  velocity  between  stations  K and  I,  The  uncertainty  in  shock 
Mach  number  from  this  calculation  should  be  less  than  1)(. 

Pres  sure  Measurements 


As  an  additional  check  on  the  shock  tube  performance  in  each  experi- 
ment, two  Kistler  type  601-B,  piezoelectric  pressure  transducers  (accelera- 
tion-compensated) are  used  to  monitor  the  pressure  variation  with  time  in  the 
driver  and  driven  sections  respectively.  The  transducer  output  in  each  case 
is  taken  directly  to  a charge  amplifier  (Kistler,  Model  504)  with  the  final 
signal  displayed  on  a designated  oscilloscope  as  shown  in  Fig.  3. 

The  combustion  dblVer.pressure  transducer  is  mounted  on  the  closed 
end  of  the  driver  section  end  given  a 2 mm  thick  surface  coating  of  silicone 
rubber  (G.E.  RTV)  for  additional  protection  from  the  higher  sustained  com- 
bustion temperatures,  A typical  combustion  trace  is  shown  in  Fig.  5a  in 
which  the  sweep  is  triggered  as  the  capacitive  discharge  began  through  the 
tungsten  ignition  wire.  The  first  portion  of  the  trace  remains  level  at 
the  pre -combustion  mixture  pressure  (423  psi  in  this  case)  where  the  original 
base-line  position  before  introduction  of  driver  gas  has  been  indicated. 

(The  measurement  of  mixture  pressure  is  only  possible  by  using  the  long  time 
constant  setting  of  the  charge  amplifier  and  can  be  maintained  for  only  a 
period  of  perhaps  3 minutes).  In  this  experiment,  a time  of  about  4 milli- 
seconds elapsed  before  the  wire  was  sufficiently  hot  to  initiate  ccmbustion. 
The  small  perturbation  of  the  signal  at  the  beginning  of  the  pressure  trace 
is  believed  to  be  noise  as  a result  of  the  electrical  discharge  throuj^  the 
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wire  and  not  a real  pressure  pulse,  particularly  since  several,  perturbations 
of  this  kind  were  observed  whenever  the  wire  was  found  to  have  broken  during 
combustion. 


The  pressure  is  then  seen  to  rise  sBDOothly  until  the  arrival  of 
the  nonstationary  rarefaction  wave  indicating  the  breakage  of  the  diaphragm. 

In  this  e3q)eriiiient,  the  diaelnragm  was  scribed  to  burst  at  a pressure  of 
24^2  psi,  while  the  combustion  pressure  at  the  closed  end  appeared  to  reada 
about  2500  psi,  possibly  indicating  a sli^t  "over-driving”  of  2^.  This  is 
probably  a result  of  the  10^  safety  factor  used  to  ensure  diaphragm  burst 
as  discussed  previously,  and  appears  inconsequential  to  the  desired  results. 

The  conbustion  pressure  measurements  are  used  to  ensure  that  correct  driver 
performance  is  obtained  in  each  e3q>eriment,  and  results  substantiated  that 
by  following  the  careful  driving  techniques  outlined  in  Section  2.1, 
excellent  control  and  repeatability  could  be  obtained. 

The  driven  section  pressure  transducer  is  flush-mounted  in  the 
selected  instrumentation  port  (Station  M in  Fig.  3) • Because  of  the  extremely 
rapid  pressure  change  across  the  shock  front,  a significant  amount  of  noise 
is  superliqposed  on  the  dutput  signal  by  the  oscillation  or  "ringing"  of  the 
transducer  diaphragm  whose  resonant  frequency  t^peared  to  be  around  1^0  Kc, 
consequently,  the  ndise  is  reduced  by  passing  the  amplified  signal  through 
a 100  Kc  filter,  although  this  limits  the  rise-time  to  about  10  microseconds 
and  hence  excludes  measurements  near  the  shock  front  itself. 

Figure  5b  shows  the  pressure  trace  obtained  for  a shock  wave 
travelling  at  Mach  nuaber  of  15  into  an  initial  argon  pressure  of  about  6 
torr.  The  small  oscillations  on  the  trace  are  the  filtered  result  of  the 
aforementioned  "ringing"  and  are  more  Intense  for  stronger  shock  waves  as 
expected.  Using  the  velocity  measurements  to  determine  shock  Mach  number  at 
station  M,  the  frozen  (liaaediately  behind  the  shock  front)  and  equilibrium 
pressures  predicted  from  the  Hugonlot  shock  Jump  conditions  are  indlceted. 

It  is  seen  that  the  pressure  soon  after  the  initial  rise  is  predicted  well 
by  the  equilibrium  shock  Jump  conditions. 

The  Initial  frozen  pressure  is  not  measured  accurately  in  comparison, 
because  of  the  flMte  rise- time  of  10  iisec.  However,  the  fact  that  a 
relaxation  time  does  exist  in  which  the  pressure  changes  from  a frozen  to 
an  equilibrium  value  is  seen  by  a change  in  the  slope  of  the  initial  pressure 
rise.  In  Fig.  5b  a small  "kink"  is  visible  as  the  pressure  appears  to  rise 
over  a time  of  approximately  10  usee.  Accurate  Interferometric  measurements 
made  simultaneously  showed  the  actual  reiixation  time  to  be  only  5 Msec  for 
this  experiment  (Mg  » 15)  which  is  less  than  the  rise- time  of  the  system 
and  hence  cannot  be  resolved. 

The  pressure  is  seen  to  remain  essentially  constant.  A decrease 
in  presswe  would  be  an  indication  of  the  arrival  of  the  contact  surface, 
or  the  reflected  rarefaction  wave  from  the  driver. 

The  accuracy  of  the  driven  pressure  measurements,  allowing  for 
gauge  calibration  factor,  charge  leakage  and  ringing  is  estimated  to  be 
arovmd  lOfa,  The  traces  in  each  experiment  are  used  principally  to  verify 
the  shock  jump  conditions  calculated  fremm  the  measured  shock  Mach  number, 
indicate  the  available  test  time  (and  possibly  the  relaxation  time),  and 
generally  ensure  thAt  conditions  through  the  shock  wave  are  as  expected. 
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2.3  Aceura<y  in  Measurements 


In  general,  measurements  are  divided  to  direct  and  indirect.  The 
temperature , for  exaa^le , is  measured  with  a standard  mercury  bulb  thermo- 
meter and  hence  it  is  being  measured  directly,  while  the  shock-wave  velocity 
is  obtained  through  time  and  distance  measurements,  and  therefore  an  indirect 
measurement  method  is  applied. 

While  Indirect  measurements  the  errors  arise:  strictly  from  the  resolu- 
tion of  the  measuring  device,  and  hence,  are  very  easy  to  estimate  (usually 
given) , the  approximation  of  the  error  involved  in  an  indirect  measurement 
is  more  oongplieated.  In  the  following  the  absolute  errors  involved  in  measuring 
the  incident  shock -wave  Mach  n\anber  M and  the  initial  pressure  P , are 
evaluated.  Note  that  these  errors  are  the  maximum  possible  errors,  and  hence 
in  general  the  actual  errors  are  ffiuah..lbsa;: . 

Incident  Shock-Wave  Mach  Number 

The  incident  shock-wave  Mach  number  is  defined  as 

M - ^ (2.1) 

8 It 

where  V is  the  shock-wave  velocity  and  a is  the  speed  of  sound  ahead  of  it. 
Thus,  tfie  accuracy  in  calculating  Mg  depends  on  the  accuracy  of  measuring 
V.  and  a.  However,  V and  a cannot  be  measured  directly,  they  are  obtained 

from 


V - J-  (2.2) 

8 w 

a - (2.3) 


where  s is  the  distance  between  two  shock  detectors,  t is  the  measured  time 
in  vhich  the  shock  wave  travels  the  distance  s,  and  T is  the  temperature  of 
the  quiescent  gas  ahead  of  the  shock  wve. 

The  uncertainty  associated  with  the  direct  measuring  of  s,  t<>and 

T are: 

dT  « ♦ 0.1°  K 
-6 

dt  » + 1 X 10  sec. 


dS  ■ + d/2  ■ + 0.3  cm 


where  d is  the  diameter  of  the  shock  detectors. 

Consequently,  the  absdlute  errors associated  with  these  measurements 

are: 
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E(T0  » 0.2  K 

-6 

E(t)  = 2 X 10  sec 
E(s)  s 0.6  cm 

where  E(^)  Indicates  the  absolute  error  associated  with  neasuresaent 

Using  the  elementary  error  approximations,  one  obtains  from  Egs. 
2.1  to  2.3. 


V 

E(M^)  - 7E(V  ) + -|•E(a) 
s a s d 

cl 

(2.U) 

E(V„)  » iE{8)  + 5-E{t) 

St  ^2 

(2.5) 

SisL  ■ i_  j;(t) 

a 2T  ^ ' 

(2.6) 

Inserting  Eqs.  2,3  and  2.6  into  Eq,.  2.^  yields 


. iiii.Slii.  1 Sffi. 

M.  f t 2 T 


(2.7) 


Equation  2,7  can  be  rewritten  in  a slightly  different  form 


E(Mg) 


^ + Mg  J ®(t)  + J ^ 


(2.0) 


IVplcal  values  of  T and  a are 


T • 300K 

a ■ 321  m/sec  for  argon 
a * 223  m/sec  for  krypton 
a • 326  a/sec  for  oxygen 
a ■ 351  m/sec  for  nitrogen 

Taking  the  shortest  8 for  which  the  shock  tube  velocity  is  measured,  to  get 
the  largest  possible  error,  and  insorting  the  appropriate  values  into  Eq. 
2.8  result,  in: 

E(M  ) 

■ ■ (1.05  Mg  + 10,10)  X lO”'^  for  argon 

8 
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E(M  ) 

(0.73  + 10.18)  X 10“ 

n s 

8 

E(M  ) 

-TCr-  = (1.08  + 10.18)  X 10“ 

M 8 

8 

E(M  ) 

— XT-"  (1-15  + 10.18)  X 10“ 

8 


for  krypton 


for  oxygen 


for  nitrogen 


Note  that  the  relative  error  E(M  )/M  increases  as  M increases. 

B 8 D 

Initial  Pressure  P 

' o— 

Initial  pressures  in  the  range  0.2  < P < UO  tor.  are  measured  by 
an  oil  manometer.  The  pressure  is  obtained^from  the  rei.ation 

p ..  X H(mm) 

^0  " 13 .'5951 

where  H is  the  oil-hdi'ght  difference  measured  in  the  manometer  p ..  is  the 
density  of  the  oil  and  it  is  given  by: 

P^^^  ■ 1.0690  + 9.5  X 10’^  (25-T)  (2.10) 

T is  the  oil  temperature  in  centigrad.  The  oil  temperature  is  not  measured 
directly,  however,  it  is  assumed  that  it  is  e(;ual  to  the  room  temperature 
near  it,  (T  ■ Tj^), 

Using  the  elementary  error  approximation,  Eqt*.  2.9  and  2.10  can  be 
rewritten  as 


B(P^)  • 

1 

|^H(mm)E(p^j^)  +'  E[H(mm)l  j- 

(2.11) 

- 9.5  X 

10"**  E(Tj^) 

(2.12) 

Inserting  Eq.  2,12  into  Eq.  2.11  results  in: 

E(Po)  • 133551  {9.5  X 10“**  H(mm)  E(Tj^)  ♦ P^^^  E him)]  j-  (2.13) 
Substituting  the  following  values. 
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E(Tj^)  » 0,2°K 

= 1.069  g/ca^ 

E(H)  3:  1 om 

Into  Eg,.  2.13  results  in: 

E(P^)  • 1.398  X 10“^  H(inm)  + T.863  x 10"^  (2.l4) 


The  range  of  pressures  under  consideration  (0.2  <P  < ^0  torr)  corresponds 
approximately  to  2 < H ^ 500  mm»  thus  the  maximum  errors  associated  vith 
initial  pressures  o?  0,2  and  ItO  torr  are  O.O79  torr  and  O.O86  torr, 
respectively.  For  a laboratory  tenqperature  of  2^  Egs.  2.9  and  2.10  can  he 
used  to  reduce  Eg.  2,lh  to: 

E(P^)  - 1.7.78  X 10“*^  P + 7.863  X 10”^  (2.15) 

0 0 


or 


+1.778x10“^  (2.16) 

0 o 

It  is  seen  from  these  equations  that  vhllethhe  absolute  error  increases 
as  P^  increases  (Eq.  2.1^)  the  relative  error  decreases  (Eq.  2,16). 

2,k  Mach-Zehnder  Interferometer 


The  principal  diagnostic  tool  used  with  the  UTIAS  Hypervelocity 
Shook  Tube  is  a 2S  cm  (8”)  field  of  view  Mach-Zehnder  Interferometer. 
Details  of  the  structure  and  design  of  this  particular  instrument  were 
given  by  Hall  (Hef.  8),  while  many  excellent  reports  have  detailed  the 
general  theory  and  operation.  A brief  description  of  the  apparatus  will 
be  presented  here  to  assist  in  discussions  of  spatial  resolution  and  aligns 
ment . 

General  Optics  , 

The  general  configuration  and  principal  dimensions  of  theilnter- 
feroffloter  are  shown  in  Fig.  6.  The  interference  optics  consist  of  two 
beam  splitters  and  two  full-reflecting  first  surface  mirrors,  each  23  cm 
in  diameter,  arranged  to  give  a 30°  angle  of  incidence  vith  respect  to  the 
central  li|^t  ray.  To  control  spacing,  orientation  and  focussing  of 
fringes,  remote  operation  is  provided  for  the  rotation  of  mirror  and 
splitter  82  about  two  perpendicular  axes  in  their  reflecting  surfaces. 

The  translitlon  of  splitter  S.  in  a direction  parallel  with  the  incoming 
light  beam  is  controlled  in  a^similar  manner  to  enable  path-length  matching 
between  the  two  arms  of  the  interferometer  about  the  central  order. 


-13- 


The  light  source  is  collimated  into  parallel  beams  by  a 23  cm 
diameter,  152.4  cm  (6oV)  focal-length  parabolic  mirror  P^.  A similar 
parabolic  mirror  Pp  serves  to  refocus  the  emergent  light  and  form  a part 
oJ  the  camera  system  in  focussing  the  object  plane  and  interference 
fringes.  Two  small  plane  mirrors  actually  direct  the  light  in  and  out 
of  the  interferometer  housing. 

The  interferometer  is  placed  around  the  shock  tube  with  the  test 
section  located  in  the  upper  arm,  and  positioned  to  provide  a Kinder 
arrangement,  in  which  the  distances  a and  b shown  in  Fig.  6 are  the  same. 

In  this  way,  fringe  spacing  and  orientation  may  be  controlled  mainly  by 
mircor  K-  with  little  effect  on  ;<rlnge  focussing.  A matching  cylindrical 
compensating  chamber  is  located  in  the  lower  arm,  enclosed  by  two  windows 
matching  those  in  the  test  section  in  both  dimension  and  tolerance.  This 
chasiber  Is  connected  to  a vacuum  port  in  the  test  section,  30  cm  downstream 
of  the  Interferometric  windows,  and  is  generally  maintained  at  identical 
conditions  with  the  test  section  prior  to  each  experiment.  This  is  not 
necessary,  however,  in  view  of  the  long  coherence  length  of  the  laser 
light  source. 

Pulsed  Laser  Light  Source 

A moat  suitable  interferometric  light  source,  particularly  for 
plasma  studies,  is  a TRO,  Model  104,  pulsed  laser  system,  capable  of 
operating  with  lasing  elements  (rods)  of  either  ruby  or  neodymium-doped 
glass  by  a correct  choice  of  end  reflectors.  Ihe  ruby  laser  operates  at 
a fundamental  wavelength  of  6943  A with  a spectral  bandwidth  of  less  then 
0.1  A and  the  neodymium  laser  at  10600  A with  a bandwidth  of  about  50  A. 

In  ."normil-mode"  operation,  the  lasing  element  is  pumped  by  an  Intense 
burst  of  light  from  a xenon  flash  lamp  lasting  approximately  1.2  milli- 
seconds, during  which  time,  a aeries  of  many  irregular  laser  pulses  are 
emitted  (providing  the  criteria  for  lasing  are  satisfied) . 

To  obtain  one,  large  short-duration  pulse,  the  laser  is  Q-swltched 
using  a Pockels  cell  (Isomet,  type  415A}  electro-optic  shutter  and  polarlxer, 
suitably  aligned  in  the  laser  cavity.  The  shutter  is  ''opened”  by  a Baird 
Atomic,  JM-IA  high-voltage  pulser  which  effectively  removes  the  quarter-wave 
retarding  volt age » applied  to  the  Pockels  cell  (in  this  case,  3 Kv  for  6943  A), 
in  a time  less  than  10  nanoseconds.  In  the  pulser  unit,  a variable  delay 
between  trigger  input  and  high  voltage  removal  assists  in  accurate  time 
control,  with  a synchronous  monitor  pulse  emitted  when  the  voltage  is  actually 
removed.  As  will  be  discussed,  a disadvantage  of  this  unit  is  the  fact  that 
the  voltage  is  re-qpplled  to  the  Pockels  ce-lli^tBra  lelatlvely  long  time 
of  1 millisecond,  over  which  the  shutter  remains  effectively  open  to  some 
extent.  For  faster  switching  of  high  voltage,  however,  more  sophisticated 
electronics  would  be  necessary. 

The  Pockels  cell  shutter  inhibits  lasing  during  the  flash-lamp 
pumping  cycle  until  the  ultimate  population  inversion  is  obtained  in  the 
laser  element,  at  which  time,  ideally,  the  shutter  should  be  opened  to 
produce  the  sx)st  intense  pulse.  Since  the  flash-lan^  radiation  pulse  is  not 
ideally  square,  the  best  time  appeared  to  be  around  oOO  microseconds  into 
the  cycle.  However,  with  the  shutter  still  open,  subsequent  normal-mode 
lasing  can  occur  if  sufficient  energy  is  still  available  from  the  flash-lamp. 
Consequently,  an  optimum  time  for  Q-switching  is  required  in  which  the  large 
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pulse  is  Sufficiently  Intense  to  expose  the  photographic  plates  but  far 
enough  along  in  the  puoplng  cycle  to  prevent  any  significant  "post-lasing” 
from  further  exposing  the  plates. 

For  the  ruby  laser,  a time  of  900  microseconds  was  found  experi- 
mentally to  be  appropriate.  Test  measurements  indicated,  for  this  time, 
an  output  pulse  width  of  30  nanoseconds  with  a total  energy  of  0«8  Joules 
using  a photodiode  and  thermopile trespectively,  inplying  a typical  power 
of  2k  megawatts.  However,  for  the  neodymium  laser,  a suitable  time  was  not 
found  for  which  the  intensity  was  sufficient  to  adeq.uately  expose  magnified 
Interferograms  while  preventing  a post-lasing  blur  in  the  flow  case. 
Consequently,  for  the  two -wavelength  Interferometry  the  ruby  laser  is 
adopted. 


^e  laser  is  equipped  with  a critically  oriented,  KDP  (potassium 
dl-hydrogsn  phosphate)  crystal,  mounted  directly  at  the  output  end  of  the 
laser  cavity  for  second  harmonic  power  generation.  The  careful  alignment 
of  this  ozystal  with  the  laser  is  discussed  by  Bristow  (Ref.  2),  and  the 
theory  for  second**harmonic  generation  is  discussed  in  the.  literature . 

Power  generation  for  this  process  is  typically  ? to  lOjt  efficient  for  the 
3 cm  long  crystal,  and  produces  coherent  radiation  at  3^71*5  a (with  a 
bandwidth  less  than  0.1  a)  which  is  both  parallel  to  and  concurrent  with 
the  incident  laser  radiation.  Consequ6n'*’'ly , by  Q-swltchlng  to  obtain 
sufficient  power,  the  pulsed  laser  becomes  a coherent  source  at  two  discrete 
wavelengths,  significantly  different  and  ideal  for  the  method  of  two- 
wavelength  interferometry. 

The  pulsed  laser  beam  is  focussed  with  a suitable  condensing  lens 
onto  a circular  iris,  0.8  mm  in  diameter,  located  exactly  at  the  focal  plane 
of  the  collimating  parabolic  mirror  P^  of  the  interferometer.  In  this 
manner,  the  ratio  of  light  source  diaSeter  (0.98  cm)  to  collimating  element 
focal  length  (15£<^  cm)  for  the  intc^rferometrio  system  is  kept  below  5.3  k 
10”^,  corresponding  to  a beam  divergence  angle  of  0.03°  (~2').  Since  the 
distance  over  which  the  light  travels  through  the  test  section  is  10,l6  cm, 
a limiting  spatial  resolution  of  the  order  of  0.05  mm  is  indicated.  (Note 
that  this  is  well  Within  the  boundary  layer  thicknesses  which  are  typically 
of  the  order  of  2 mm) . 

The  condensing  lens  is  chosen  from  simple  geometry  to  concentrate 
the  light  source  energy  over  the  specific  region  of  Interest  in  any 
experiment  as  shown  in  Fig.  7.  This  was  particularly  isportant  in  maintaining 
a sufficient  Intensity  for  interferograms  in  which  the  camera  magnification 
is  increased.  As  shown,  a suitable  focal  length  f , for  the  condensing 
lens  may  be  determined  from  ° 


where  d is  the  diameter  of  the  pulsed  laser  beam,  found  to  be  8 mm,D  is  the 
diameter  of  the  final  collimated  light,  and  f *>  152. U cm  is  the  focal  length 
of  the  collimating  parabolic  mirror.  ' 

Consequently,  a 22  cm  focal-length  condensing  lens  is  adeqxiate  in 
illuminating  a 6 cm  diameter  region  for  the  boundary  layer  study,  while  a 
5.3  cm  focal -length  lens  is  used  when  examining  the  entire  20  cm  diameter 
field  of  view. 


-15- 


It  Is  pointed  out  here  that  a condensing  lens  Is  always  necessary 
with  the  high  Intensity  laser  In  preventing  radiation  damage  to  the 
expensive  optical  components  of  the  Interferometer.  The  small  input  plane 
mirror  to  the  interferometer  was,  in  fact,  burned  occasionally  because  of 
the  small  beam  diameter  at  that  point,  particularly  with  longer  focal-length 
condensing  lenses. 

The  advantages  of  the  pulsed  laser  as  a light  source  may  be  summarized 
as  follows: 

(a)  The  short,  30  nanosecond  light  pulse  is  sufficiently  fast  to 
freeze  the  shock  motion  and  associated  phenomena. 

(b)  The  Intensity  of  light  at  discrete  wavelengths  over  this  short 
duration  is  sufficient  to  expose  standard  photographic  plates, 
and  is  much  hiither  than  the  background  plasma  radiation. 

(c)  The  two  discrete,  significantly  different  wavelengths  of  69^3  £ 
and  3^71  a are  suitable  for  the  resolution  of  species  densities 
in  a plasma  using  two-wavelength  Interferometry. 

(d)  The  highly  coherent,  monochromatic  beam  ensures  a long  coherence 
length  such  that  fringes  of  good  contrast  are  visible  even  for 
relatively  large  mls-match  in  optical  path  lengths  between  the 
test  section  and  compensating  arms  of  the  interferometer. 

(e)  The  beam  is  highly  unidirectional  and  can  be  easily  focussed 
down  to  a small  size  to  provide  excellent  spatial  resolution 
for  the  bounary  layer  studies. 

Light  Source  Timing 

Some  prior:  knowledge  of  the  expected  shock  speed  (at  least  within 
10J()  is  required  to  synchronize  the  pulsed  laser  operation  with  the  shock* 
tube  flow.  This  is  necessary  not  only  in  controlling  the  time  at  which 
the  laser  is  Q-swltched  to  expose  the  interferogram  plates,  but  also,  prior 
to  this , in  flash-lamp  pumping  the  lasing  element  for  the  optimum  length  of 
time.  The  calculations  required  are  most  easily  visuaxized  on  an  x-t 
diagram  such  as  Fig.  4,  With  the  time  or  shock  position  specified  at  which 
the  interferogram  is  to  be  taken,  the  closest  shock  detector  upstream  is 
designated  to  trigger  the  Pockels  cell,  with  the  time  difference  accounted 
for  in  the  high  voltage  pulser  delay  unit.  Similarly,  the  closest  shock 
detector  at  least  900  microseconds  \ipstream  is  chosen  to  Initiate  the  flash- 
lamp,  with  the  excess  time  greater  than  900  (xsec  accounted  for  in  another 
delay  unit. 

Thib  control  of  laser  light  source  operation  is  included  in  Fig.  3. 
The  laser  flash-lamp  capacitor  bank  is  usually  triggered  by  the  shock  arrival 
at  station  F,  suitably  delayed  in  a Tektronix,  ^I^pe  Oscilloscope.  The- 
Pockels  cell  Q-swltch  is  always  triggered  from  station  I,  delayed  correctly 
in  the  pulser  delay  unit,  to  take  the  interferograms  at  the  desired  time. 

As  mentlcmed,  the  time  at  which  the  interferograms  are  taken  is  always 
recorded . To  monitor  the  actual  flaah-lamp  pumping  time  in  each  experiment , 
synchronization  pulses  from  both  the  flash-lamp  capacitor  bank  and  the 
Pockels  ceil  pulser  unit  are  used  to  start  and  stop  respectively  a UTIAS 
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microBecdnd  counter.  It  was  found  that  variations  of  ^ to  100  ijpec  in  this 
time  did  not  seriously  affect  either  the  quality  of  intensity  of  the  inter- 
ferograms,  thereby  allowing  for  some  error  in  the  actual  shock  speed  from 
that  predicted,  for  the  long  distance  involved. 


Camera  System. 

The  test-section  object  plane  that  is  defined  by  cross-wires 
installed  in  the  test  model  or  on  both  windows,  is  focussed  onto  fi^  plates 
or  ground  glass  screens  by  the  emergent  parabolic  mirror  Pg  of  the  interfero- 
meter followed  by  one  or  two  converging,  achromatic  lenses,  the  number  depending 
on  the  desired  magaification.  The  parabolic  mirror  is  located  approximately 
206  cm  from  the  object  plane  and  focusses  the  image  of  the  plane  at  a distance 
of  550  cm  from  the  mirror  with  a corresponding  magnification  of  2.05.  To 
record  interferograma  covering  the  entire  20  cm ’.'field  of  view,  an  over^l  ^ 
magnification  of  1/2  is  necessary  when  using  standard  10  cm  x 13  cm  (ft’  x 5 ) 
film  plates.  As  shown  in  Pig.  7«  this  is  accomplished  by  inserting  an  01  cm 
focal-length  lens,  for  convenience,  near  the  focal  plane  of  the  parabolic 
mirror* 


To  obtain  an  Image  magnification  of  3 (needed  for  boundary  layer 
studies  for  example),  the  camera  assembly  is  merely  moved  back  and  an  additional 
10  cm  fbcal-length  lens  inserted,  approximately  21  cm  past  the  focussed  image 
from  the  first  lens  as  shown  in  Pig.  7b.  The  actual  distance  between  lenses 
Is  sensitively  adjusted  in  focussing  the  camera  and  setting  the  overall 
maffiificatlon.  It  should  be  pointed  out  that  the:  lenses  were  selected  to  be 
physically  compatible  with  the  camera  assembly  and  to  carefully  avoid  laser 
beam  focussing  on  any  camera  components. 

In  the  camera  assembly  itself,  simultaneous  interferograma  at  two 
wavelengths  are  obtained  using  a small  beam  splitter  and  plane  folding  mirror, 
followed  by  respective  line  filters  to  separate  the  con^ponent  wavelengths 
and  shield  the  plates  from  plasma  radiation.  Kodak,  Boyal  X Pan  (1250  ASA), 

10  cm  X 13  cm  plate  film  are  used  to  recoi-d  the  interferograms.  Plates  are 
subsequently  developed  (7  minutes)  in  Kodak,  DK  50  solution  for  best  results. 

The  mechanical  shutter  of  the  camera  is  opened  immediately  before  recording, 
allowing  the  30  natiosecond  laser  pulse  to  expose  the  plates.  This  is  necessary, 
of  course,  since  no  mechanical  shutter  coxild  react  to  flow  changes,  typically 
of  the  order  of  5 mm/(tfiec. 

Critical  Alignment  Procedures 

In  the  two-dimensional  evaluation  of  interferograms  it  is  generally 
required  that  the  collimated  light  of  the  interferometer  test  beam  be 
perpendicular  to  the  plane  in  which  flow  changes  occur,  or,  in  effect,  parallel 
with  the  surface  over  which  the  flow  travels.  !lhis  condition  of  no  changes  in 
the  lie^t  path  direction  enables  the  integrated  effect  of  the  optical  path 
length  on  the  interference  to  be  resolved  directly  into  its  components  of 
refractive  index  and  distance.  Thus,  any  misalignment  constitutes  an  error  in 
the  resulting  measurement.  This  becomes  particularly  iitportant  in  boundary  layer 
studies  where  significant  flow  changes  occur  within  a region  perhaps  2 mm  from 
the  flow  surface, and  a high  degree  of  spatial  resolution  is  necessary.  In  this 
case,  misalignment  not  only  causes  a severe  loss  of  accuracy  but  often 
prevents  measurexaents  otherwise  possible  near  the  surface  because  of  surface 
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reflection  and  diffraction.  Consequently,  special  care  is  taken  to  sensitively 
align  the  interferometer  with  the  teat  model  surface  and  optimize  the  spatial 
resolution . 

As  mentioned,  the  effective  size  of  the  laser  light  source  was  made 
small  to  maintain  good  spatial  coherence  across  the  field  of  view  and  thus 
ensure  that  the  beam  is  sufficiently  parallel  within  itself  to  permit  a 
resolution  of  about  0.0^  mm.  Alignment  of  this  test  beam  with  the  test  model 
surface  is  then  carried  out  using  a procedure  which  ultimately  made  use  of  the 
spufious  interference  fringes  which  form  adjacent  to  an  extended  solid  surface 
when  slight  misalignment  is  present.  This  method  was  first  proposed  by  Howes 
and  Buohele  (Ref.  9)t  and  has  since  been  used  in  other  studies. 

In  view  of  the  large  physical "size  of  the  interferometer  (Fig.  6), 
initial  steps  are  taken  in  order  to  facilitate  the  alignment  procedure  whenever 
necessary.  A platform  was  built  and  attached  to  the  housing  of  the  Interferometer 
on  vhlch  the  appropriate  light  source  could  be  positioned.  In  this  manner., 
the  light  source  and  Interferometer  move.-i  as  a unit  when  alignment  is  carried 
out.  Tne  focal  plane  of  the  collimating  parabolic  mirror  P.  'is  then  determined 
by  first  attaching  a small  light  bulb  behind  a circular  Irll,  adjustable  in 
diameter,  and  directing  this  light  into  the  interferometer.  By  aligning  the 
input  plane  mirror  and  the  parabolic  mirror  P.,  and  moving  the  iris  (and  light 
bulb),  a position  vis  found  over  the  platform  where  the  reflection  of  the  light, 
returned  from  the  windows  of  the  compensating  chamber,  focussed  directly  back 
to  the  centre  of  the  iris . The  iris  Itself  vis  rigidly  affixed  to  the  light 
source  platform  at  this  position,  effectively  defining  a focal  point  of  the 
parabolic  mirror  P. . In  this  manner,  it  is  ensured  that  any  light  directed 
through  the  iris  would,  on  reaching  the  parabolic  mirror,  form  a parallel  beam 
which  would  be  perj^ndidular  to  the  eompensating  chamber.  It  should  be  pointed 
out  that  this  iris  position  was  never  altered  again. 

To  align  the  interferometer  with  the  shock  tube  bottom  or  the  surface 
of  the  test  model,  a small  heliw-neon  laser  (Spectra  Physios,  Model  13^) 
is.i  mounted  on  the  platform.  By  directing  the  laser  beam  through  the  0.8  mm 
diameter  iris  from  various  angles,  it  ^is  verified  that  the  reflection  of  the 
beam  from  the  compensating  chamber  windows  always  returned  to  the  centre  of 
the  iris.  The  entire  interferometer  as  a unit,  including  cooqpensatlng  chamber 
and  light  source,  is  then  adjusted  in  both  pitch  and  yaw  with  the  test  section 
of  the  shock  tube  until  the  leflections  of  the  helium-neon  laser  beam  from 
the  test  section  windows  are  also  returned  to  the  centre  of  the  iris.  !Ihis 
provides  a good  approximation  to  the  correct  alignment  with  the  shock  tube 
bott(»i  or  with  the  surface  of  test  models  since  the  test  section  windows  wave 
very  nearly  perpendicular  to  the  shock  tube  surfaces.  A more  sensitive 
procedure  is  then  carried  out. 

To  illuminate  the  teit  section,  the  2 mm  diameter  gas  laser  beam  is 
expanded  by  inserting  a short,  1 cm  focal  length  "eyepiece"  lens  in  the  same 
maxmer  as  shown  in  Fig.  7,  eu^  that  the  beam  initially  focussed  at  the  iris. 

A moveable  ground-glass  screen  is  placed  at  the  image  plane,  formed  by  the 
parabolic  mirror  Pg,  of  the  cross-wires  d'^coussed  previously.  The  inter- 
ferometer is  then  adjusted  to  focus  many  vertical  fringes  on  the  screen,  in 
a desired  direction.  (It  is  pointed  out  that  fringe  spacing,  etc,  are 
controlled  by  mirror  Mg  and  splitter  Sg  of  the  interferometer,  which  do  not 
affect  the  incoming  beam  alignment  with  the  test  section.)  Using  a small 
"Jeweller's"  eyepiece  behind  the  screen,  the  interference  immediately  adjacent 


to  the  surface  can  then  be  observed  at  a magnification  of  about  20. 

As  predicted  by  Howes  and  Buchele  (Ref,  9),  spurious  diagonal  fringes 
are  observed  at  the  surface  when  the  beam  is  misaligned.  When  the  incoming 
beam  is  inclined  at  a small  negative  angle  with  respect  to  the  surface  these 
spturious  fringes  result  from  light  that  is  reflected  from  the  model  surface 
(and  also  diffracted  by  the  surface  edge  nearest  the  screen)  and  interfering 
with  the  reference  beam.  When  the  beam  is  at  a small  positive  angle, 
diagonal  fringes,  oppositely  inclined,  resulted  from  the  light,  diffracted 
by  the  surface  edge  nearest  the  laser  light  source.  Interfering  with  the 
reference  beam.  When  correct  alignment  is  obtained,  these  spurious  fringes 
effectively  disappeared,  althoi^gh  a slight  '•picket-fence'*  effect  is  observed, 
caused  by  the  combination  of  the  two  opposed  patterns  and  the  diffraction 
pattern  produced  by  the  surface  edge  nearest  the  light  source . 

Consequently,  while  observing  this  interference  pattern  on  the  screen, 
the  entire  Interferometer  is  sensitively  adjusted  in  pitch  (usually  by  air 
pressure  in  the  pneumatic  tire  mounts)  until  the  correct  pattern  is  obtained. 
The  final  alignment  is  generally  approached  from  the  negative  angle  of  the 
incoming  beam  with  the  surface,  since  the  spurious  fringes  formed  by  the  re- 
flected light  are  more  distinct  than  those  formed  by  diffraction,  mxe  correct 
alignment  with  respect  to  yaw  is  also  found  by  adjusting  fringes  perpendicular 
to  the  downstream  surface  of  the  test  model  and  using  the  spurious  fringes 
formed  adjacent  to  this  surface  as  well. 

Bunting  and  Devoto  (Ref.  10)  found  that  the  residual  misalignment 
following  this  type  of  procedure  was  less  than  the  uncertainty  of  the  beam 
collimation  arising  from  the  finite  size  of  the  light  source.  Consequently, 
the  limiting  spatial  resolution  after  alignment  should  remain  as  0.05  mm. 

It  should  be  noted  that , after  this  alignment  i«  carried  out , the 
reflections  of  tue  laser  beam,  returned  from  the  test  section  windows,  are- 
found  to  be  displaced  downward  from  the  centre  of  the  iris  by  about  U mm. 

This  indicates  from  simple  geometry  (for  the  152. U cm  focal  length  collimating 
parabolic  mirror  F^)  that  the  test  section  windows  'are'  about  5'  of  arc  from 
being  perpendiculiu'  to  the  flow  surfaces.  The  position  on  the  iris  where 
these  reflections  'are  returned  have  been  recorded  to  serve  as  an  additional 
indicator  as  to  the  positioning  for  correct  alignment. 

Following  this  procedure,  final  adjustments  are  made  to  the  system. 

The  kelium-neon  laser  is  removed  from  the  light-source  platform  and  placed 
on  a stand  behind  -the  platform.  The  narrow  laser  beam  is  again  directed 
into  t^e  Interferoneter  through  the  iris  so  as  to  pass  through  the  centre 
of  Interest  in  the  test  section.  The  camera  system,  including  lenses  and 
film-plate  holders  are  then  positioned  along  the  axis  of  '::his  beam. 

Similarly,  the  pulsed  laser  light  source  is  then  mounted  on  the  platform 
between  the  gas  laser  and  the  iris  and  adjusted  axially  with  the  gas  laser 
beam.  (This  is  easily  accomplished  by  opening  the  back  of  the  pulsed  laser, 
allowing  the  beam  to  pass  through  the  laser  element,  and  aligning  by  means 
of  cross-hairs  at  each  end  of  the  laser  cavity) . Final  camera  adjustment 
and  focussing  of  fringes,  etc,  is  done  using  a tungsten  light  placed  at  the 
iris  (with  the  iris  wide  open)  to  more  easily  define  the  correct  plane  of 
focus.  (The  white  light  source  is  also  used  to  path-length  match  the  two 
amui  of  fnt  interferometer  to  give  the  best  fringe  contrast  about  the  central 
order . ; 
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The  procedures  discussed  In  this  section  can  he  performed  easily  and 
efficiently  vhenever  necessary,  and  present;  a marked  improvement  over  methods 
previously  used  on  this  facility. 

3.  PRINCIPLES  OF  TWO-WAVELEKGTH  INTERgSROMETRY 


The  Madi-Zehnder  Interferometer  provides  a means  of  spatially  resolving 
the  phase  refractive  index  throughout  a two-dimensional  flow  field  in  the  test 
section.  In  general,  the  change  in  refractive  index  from  some  reference 
condition  may  he  detenfiined  hy  measuring  tibe  change  in  the  Interference 
pattdtAoor  .’'fringe  shift”  from  this  reference.  Specific  evaluation  techniques 
are  discussed  in  detail  in  Section  U.  Using  the  reference -point  method  of 
analyzing  Interferograma  as  descrlhed,  the  change  in  interference  is  related 
to  the  change  in  refractive  index  (Eq.  4.13) 

n(x,y)  - n(x^,y^)  « ^S(x,y)  (3.1) 

where  S(x,y)  is  the  nondimensional  "fringe  shift”  or  change  in  interference 
measured  at  the  position  (x,y)  relative  to  the  interference  change  at  a 
reference  point  (x^,y^),  as -definedot  in  Section  4,  Eq,  4.12, 

^ is  the  wavelength  of  incident  light  used, 

L is  the  width  of  the  test  section,  across  which  the  Itgt.t 
beam  travels  (I0.l6  cm  ■ 10. l6  x 10°  A°  in  this  study) 

n(x,y)  is  the  refractive  index  at  the  position  (x,y)  in  the  flow 

field,  and  n(x  ,y  ) is  the  refractive  index  at  the  reference 
point  (x^,y^)  In  the  flow,  and  must  be  known  if  an  absolute 
determination  of  n(x,y)  is  to  be  made. 

For  dilute  atomic  or  molecular  gases,  the  phase  refractive  index 
can  ge  nerally  be  related  directly  to  the  density  by  the  classical  Oladstone- 
Pale  equation, 

n-1  » KP  (3.2)' 

where 

p ■ mn  is  the  gati  density,  which  can  be  written  in  terms  of 
particle  mass  m,  and  number  density  n,  and 

K is  the  Qlads tone -Dale  constant , which  is  constant  for 
a given  wavelength  of  light  and  over  a large  range'  Pf 
pressures. 

The  Oladstone-Dale  constant  may  be  found  frdm  the  Induced  electric  dipole 
polarizability  of  the  gas,^  , as  (Ref.  2): 

K » 2Tr  (5  /m 

In  a dilute  mixture  where  the  interaction  between  constituent  gases 
may  be  neglected  relative  to  the.  interaction  with  the  electric  field 
associated  with  the  light  wave,  the  refractlvlties  are  additive  8U(^  that 


the  total  refract  Ivlty  of  the  mixture  Is  given  hy  a siumoation: 


n-1  - ^ (n-l)g  - ^ KgPg  ^3^3) 

s s 

where  the  subscript  "s"  refers  to  a species  of  the  mixture. 

Singly  Ionized  Plasma 

As  discussed  by  Alpher  and  White  (Ref.  U),  most  plasmas  may  be 
considered  in  this  manner,  at  least  for  Incident  radiation  in  the  optical 
region  of  the  spectrum.  Consequently,  for  a singly  lohlswd  plasma  consisting 
of  a mixture  of  neutral  atoms,  ions  and  electrons,  the  total  plane  refractive 
index  may  be  written  as: 


n-1  « (n-1)^  + (n-l}j^  + (n-l)^ 


(3.4) 


where  the  subscripts  a,  1,  e denote  atoms,  ions  and  electrons  respectively. 
It  should  be  note^  that  the  neutral  atom  and  ion  components  cohtaln 
Isqpllcltly  the  contributions  from  the  excited  state  as  well  as  ground  state 
species . 


The  refractlvity  of  electrons  can  be  obtained  from  the  plasma  dis- 
persion formula.  With  no  applied  magnetic  field,  and  the  frequency  for 
momentum  transfer  collisions  between  electrons  and  heavy  particles  shall 
with  respect  to  the  frequency  of  the  light,  it  may  be  shewn  that 

2 

w 

(n-l)^  « - (3.5) 

2u 

where 

n jat 

__  * 

y ■■■  ■ ■■ 

P “a 

is  the  classical  plasma  frequency,  and 
m^  is  the  mass  of  atom. 

W = 2Tr  c/X  = 1.88365  X 10^^/X 
w is  the  angular  frequency  of  the  light  source. 

Finally  in  terms  of  wavelength,  X(in  A°),  £q.  3.5  becomes: 


(n-1)^  - 


4.4849  X 10"^^ 
aa 


(3.6) 
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The  highly  dlsperolve  nature  of  the  electrons  is  easily  seen  by  the  dependence 
on  the  square  of  the  waurelength.  The  total  plasm.  i ifraotive  index  may  now 
he  written  as : 

2 

n-l  =.  + K^pj  - 11.481(9  X 10-1'*  ^ 

a 

In  a singly  ionized  plasma,  n^  = n^,  m^  total  plasma 

density  is  given  by: 


(3.7) 


Hence , rewriting  the  plasma  refractive  index  in  terms  of  p and  n : 

* 

,-30  >2 


n-*! 


. 4.4849.x  10-3. 1 ^ „ 

1 i m a e 

A •• 


It  is  convenient  to  define 

Z(X)  = ( K -K.)  + — — 
a 1 m. 


-30  ^2 


(3.8) 


The  total  plasma  refractive  index  may  then  be  expressed  as: 


n-l  « K p - Z m n 

ft  ft  ft 


(3.9) 


The  total  refractive  index  for  a singly  ionized  plasma  is  seen  to 
depend  explicitly  on  the  plasma  total  density  and  electron  number  density, 
and  inqplicltly  on  the  wavelength  of  light  through  the  known  dispersive 
constants,  K and  Z.  Consequently,  if  the  plasma  is  examined  interferometrlcally 
using  light  Sit  two  discrete  wavelengths  simultaneously,  two  values  of  the 
plasma  refractive  index  may  be  found  for  the  same  plasma  conditions.  From 
Fq.  (3.9),  two  equations  in  two  unknowns  result,  and  the  plasma  total  density 
and  electron  number  density  can  be  resolved.  This  forms  the  basis  for  the 
method  of  two-wavelength  interferometry,  whose  sensitivity  depends  on  the 
dispersive  effect  of  the  electrons.  It  is  important  to  note  that  no  assump- 
tions have  been  made  to  this  point  as  to  the  plasma  being  in  either  thermal 
or  chemical  equilibrium. 

It  is  most  common  to  relate  the  plasma  densities  directly  to  the 
interference  measurements,  using  Eq.  (3.9)  to  replace  refractive  index  in 
Eq.  (3.1): 


\ [p(x,y)-p(x^,y^)]  - Z mjn^(x,y)-njx^,y^)]  = ^ S(x,y)  (3.10) 


In  shortened  notation,  dropping  the  position  coordinates  and  the  subscript 
'*a”  from  the  Gladstone-Dtile  constant. 
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s 


K ^ - Z A Hg 


Vlheu  tvo  interferograms  are  taken  slmultaneovisly  at  two  different  wavelengths, 
two  distinct  measureaents  of  fringe  shift  are  made.  Using  subscripts  "I” 
and  ”2"  to  denote  the  two  wavelengths,  then. 


h 


Solving , 

- Wi 


W2  - Wi 


(3.U) 


(3.12) 


The  constants  in  these  eolations  for  the  two  wavelengths  of  the 
previously  discussed  rhby  laser  (X  ■ 6943  A,  Xg  ■ 3471.?  %)  for  argon  and 
kzypton  are  listed  in  Table  2.  ^ 


Note  that  the  first  term  in  E%.  3.6  (l.e.,  K -JC.)  is  at  least  one  order 
of  magnitude  smallsr  than  the  electron  dispersion  tem,^and  hence,  an  accurate 
value  for  this  ratio  is  not  critical  to  the  evaluation. 


Diatqnic  0^ 

For  a diatomic  gas  Eq.  (3.3)  can  be  written  in  the  form 

n-1' - (n-1)  - (n-l).  (3.13) 

2a  ft 


where  tt  designates  a molecule  and  a an  atom. 

Ihus  Bq.3.10  now  becomes: 

s • ^ { p [y  * VI  -Prty^-“r>  * Vr]  } 

OX  after  rearrangment 

s . £ { yiy-y  - ( W (ea-p^a,)  j-  (3.11.) 

Following  the  above  procedure  of  solving  two  equations  (two  different  wave- 
lengths) with  two  unknowns,  4p  andA(pa}  results  in; 
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(3.15) 


1 ^ Sg  - Zg  \ S; 


A(pa)  = 


h - \. 


O.iiS) 


vhere  Z^  » “'*  Gladstone-Dale  constants  for  the 

molecule  and  the  atom  res^ctlvely. 


It  is  worth  noting  that  the  final  form  for  both  singly  ionized  plasma 
(of  monatomic  gases)  and  diatomic  dissociating  gases  is  exactly  identical 
(cooqpare  E^s.  3.11,  3.12,  3.13  and  3.16).  The  constants 

and  Zg  for  oxygen  and  nitrogen  are  given  in  Table  3. 

3.1  Absolute . Errors 


%• 


The  errors  in  the  Indirect  meas\irements  of  the  total  density,  the 
nusiber  density  the  quantitative  pCX  and  (o-degree  of  dissociation,  X'^ogree 
of  ionization)  changes,  that  are  directly  related  to  the  accuracy  within  which 
the  fringe  Shifts  can  be  measured,  can  be  calcualted  using  the  elementary  error 
approximation . 

Total  Density 

Using  the  elementary  error  approximation  Eqs.  3.11  and  3.15  result  in: 


B(Ap)  = + •[  I I E(S2)  + I 11^1  } (3Ut) 


Note,  E(0.)  Indicates  the  error  associated  with  measurement  Thus  if  the 

errors  in  measuring  respective  fringe  shifts  are  about  the  same,  i.e,  E(S.)  * 
£(82)  ■ E(S),  Eq.  3.17  together  with  Eq.  3.11  or  3.15  yields 

[ 1^1^ I (3.18a) 

for  a monatomic  gas,  and 


I 1 ^2^.1 


(3.18b) 


for  a diatomic  gas. 

Using  the  constants  from  Tables  2 and  3 for  argon,  krypton,  oxygen  and 
nitrogen  respectively,  Eqs.  3.l8a  and  b reduce  to  (h'ote  L « 10. 16  cm): 


E(^ ) * 4.292  X 10  ^ E(s)  &—  t 


-24- 


E(Ap)  » 5.866  X 10“^  E(s) 

•5 

for  krypton 

aoT 

E(Ap)  « It. 031  X 10“*^  E(S) 

t- 

3 

for  oxygen  , 

cm 

E(^)  » 2.097  X 10"^  E(S) 

for  nitrogen 

cm 


A conservative  estimate  of  the  error  in  measuring  fringe  shift  Is  0.0?.  Hence 
the  estimated  errors  In  the  density  change  measurements  are: 


£(Pp) 

■ 2.146  X 10"^  g/cB? 

for  argon 

E(^  ) 

- 2.933  X 10“^  g/ov? 

for  krypton 

E(^  ) 

■ 2.016  X 10"^  g/cm^ 

for  oxygen 

E(^  ) 

- 1.049  X 10"^  g/cA^ 

for  nitrogen 

Ihe  Quantities  pOand  fiil 

Ihe  shaolute  errors  in  maasurlng  the  quantities  pa  or  0.x  can  he  written 
(from  Eqs.  3.12  end  3«l6)  in  the  form 


E[A(pA)]  •»t{\  J 


where  A stands  for  either  X or  a and  PX  " m n^ . Following  the  shove  procedure 
results  in: 


®E^(PX)  ] “ LEKgZ^-  K^Zg)  [ (3.20a) 

E[A(po)  ] *•  ]®ES)  (3.20b) 

Using  the  constants  from  Tables  2 and  3 will  give.:  ■ 

eEA(px)]  - I*. 101  X 10*^  E(S)  for  argon 

cm^ 
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far  krypton 


(PX)]  = 8.408  X 10“^  11(8) 

1 cm 

B[A  {pa)  ] = 1.229  X lO’^  E(s)  ^ for  oxygen 

E[A  (pa) 3 *•  1.221  X 10“^  E(S)  ^ — for  nitrogen 

cm? 

snd  finally,  if  E(s)  is  taken  as  0.0^ 

EC^(PX)]  = 2.051  X lO"*^  for  argon 

ClBr 

E[A(px)3  “ *+«204  X lO"*^  for  kryptcn 

CSl 

-4  <r 

E[A(pa)3  ■ 6.145  X 10  ^ for  oxygen 

ca'^ 

E[A(pa)3  “ 6.105  X lO"^  ^ for  nitrogen 


EuBfcer  Density 


be  written: 


By  the  deflxiltlon  of  nunber  density  the  following  relations  can 


n M Bit 
e a 

A 


« . i 

. f 


n - £2 
a a. 


Thus,  the  absolute  error  In  the  nunber  densities  can  be  lanedlately  obtcJ.ned 
from  the  req;>ectlve  errors  associated  with  the  quantities  pa  and  px 
(calculated  bbove) . Consequently,  the  absolute  errors  are; 


A(ng)  = 6.182  X 10^^  E(S)  -2^  for  argon 


M«-)  = 6,108  X 10^^  E(S)  for  krypton 
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A(na) 

» 4.626  X 

10^°  E(S) 

1 

3 

for  oxygen 

cm 

A(na) 

» 5.249  X 

10^^  E(S) 

1 

3 

cm 

for  nitrogen 

Finally,  if  E(S)  is  assuined  to  "be  0.05  the  absolute  errors  hesoae: 


A(n^) 
A(n  ) 
^(n  ) 

A 

A(n^) 

It  should  he 
hence,  place 


3.091  X 10^^ 

1 

3 

fo3.’  argon 

cm 

3.054  X 10^^ 

1 

for  krypton 

cm 

2.313  X 10^^ 

1 

0. 

for  oxygen 

cm*' 


• 2.625  X 10^®  ^ for  nitrogen 
cm'^ 


noted  that  all  the  above  calculated  errors  have  fixed  values  and 
a lover  limit  on  the  changes  which  can  be  measured. 


In  summary t tvo>-vavelength  interferometry  provides  an  excellent  means 
of  measuring  changes  in  the  total  density,  the  nuxaber  density  and  the  (quantity 
PCX  (or  PX  ) if  the  measured  values  are  at  3 east  one  order  of  magnitude  greater 
than  the  absolute-errors  associated  with  the  measurements,  i.e.,  in  the  case 
of  argon  p > 10 g/cm^,  pX  > 10~5  g/cm3  and  n > 10^“  cm-^.  The  measurement 
technique  does  not  perturb  the  gas  under  consideration  to  any  appreciable 
extent , and  may  be  applied  equally  veil  to  nonequilibrium  as  veil  as 
equilibrium  situations, 

1*.  SPATIAL  ANALYSIS  OF  HTTERFEROORAMS 


In  order  to  obtain  quantitative  measurements  of  gas  densities  using  an 
interferometer,  care  must  be  taken  in  reducing  the  resulting  fringe  patterns 
correctly.  As  the  flow  being  studied  becomes  more  complex,  the  simpler  direct 
evaluation  methods  benome  difficult  to  apply.  Such  an  exsjiqple  is  shown  in  the 
interferogram  (Fig  .8) of  the  shock  structure  in  argon  near  the  shock  tube-. vail. 
Interferograms  such  as  tliesQ,vlth  their  dark  and  bright  lines  of  interference, 
or  fringes,  contain  a vast  amount  of  information  about  the  refractive  index 
(and  hence  density)  variations  over  the  entire  field  of  view.  In  order  to 
extract  the  maximum  amount  of  data  in  flows  of  any  complexity,  a digital  evalu- 
ation program  was  set  up  for  use  with  the  23  cm  (9  in)  field  of  view  Mach- 
Zehnder  interferometer,  in  which  the  spatial  coordinates  (x,y)  of  the  various 
lines  of  Interference  (fringes)  on  the  Interferograms  are  put  into  digital 
form  such  that  a cooputer  analysis  is  then  possible. 

In  developing  the  relationship  between  orders  of  interference  and 
refractive  index,  the  basics  of  interferometry  are  reviewed  briefly  as  appli'e'dj' . 
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to  the  Mach-Zehnder  interferometer  vhose  con^onents  are  shovn  in  Fig.  9. 
(More  details  may  he  found  in  Refs.  12  and  13).  Light  from  a coherent 
monochromatic  source  is  divided  hy  a first  beam  splitter  S. , into  two 
coherent  beams,  one  of  ^ich  travels  through  the  test  section  p,  vhile  the 
other  travels  through  a cosg)ensation  chanber  %.  When  reunited  at  the  second 
splitter,  S2,  the  light  waves  will  interfere  constructively  or  destructively 
depending  on  the  relative  time  of  travel  through  each  arm.  The  time  taken 
for  a plane  electromagnetic  wave  to  pass  through  an  isotropic  medium  may  be 
written  in  terms  of  the  "optical  path  length”  t as 


At  ■ i-  T 

% 

where  c^  is  the  speed  of  light  in  vacuo,  and 

T ndl>^- 

vhere  n is  the  refractive  index,  which  is  a function  of  the  medium  and  of 
the  frequency  of  the  wave,  and^  is  the  geometrical  path  over  which  the 
wave  travels,  or  the  arc  length.  For  coherent  light,  constructive  inter~ 
farence  occurs  if  the  difference  in  optical  path  lengths  between  the  two 
arms  is  an  even  multiple  of  half  wavelenghs  of  the  light  used,  while  des- 
tructive interference  takes  place  if  this  difference  is  an  odd  multiple  i.e., 


r 2N’  for  maxima 

1 (2N'  + 1)  V2  for  minima  {k.i) 


where  K’  ■ 0,  + 1,  + 2,...,  may  be  termed  the  order  of  the  interference,  "p” 
denotes  a wave~traveiling  through  the  test  section  am,  "q”  denotes  a wave 
travelling  through  the  compensation  arm.  For  convenience,  let  the  optical 
path  lengths  in  the  respective  arms  be  denoted  as 


n di^ 

<1  q 


p S2  P ® p b p So 

L n d/  = L n di  + 1 n dX  + , n di 

jSi  P P jSi  p p ^ p p PI 


Here,  the  optical  path  length  through  the  test  section  has  been  broken  at 
"a",  the  entrance  position  to  the  test  section,  and  "b",  the  exit  position 
from  the  test  section,  into  three  additive  coisponents.  The  tem  t*,  is 
then  defined  as: 


T*  = 


T 


q 
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(‘♦.2) 


such  that  T*  accounts  for  all  the  optical  path  lengths  outside  the  test 
section.  The  Interference  (i^.l)  may  then  be  written  as 


r n di  + T*  = -f 

Ja  P P I 


2N'  maxima 

(2N*  + 1)  y^2.  for  minima 


To  write  an  easier  and  more  general  relationship,  a continuous  linear 
variable,  N,  Is  considered  such  that: 


n di 
P P 


+ T*  = NX 


(1^.3) 


where,  for  N ■ 0,  + 1,  + 2,...,  there  will  be  total  constructive  Interference, 
or  maxima;  while  for  N ■^  + 1/2,  ^ 3/2,  + 5/2,, , , , there  will  be  total 
destructive  Interference ,**or  minima.  For  all  other  intermediate  values  of 
N,  there  will  be  varying  "gray"  Interference,  or  a graduation  of  these  two 
extremes . 


The  quantity  of  interest  is  generally  the  refractive  Index  n of  the 
medium  in  the  test  section,  which  is  then  related  to  the  density?  Thus,  If 
the  path  teken  by  the  light  travelling  through  the  test  section  is  known,  and 
all  the  optical  path  lengths  # outside  the  test  section  are  accounted  for, 
then  orders  of  interference  N measured  from  an  interferogram  may  be  used  with 
(^.3)  to  determine  the  integrated  or  cumulative  value  of  n_  over  this  path. 
To  dete»ine  for  this  equation,  the  interference  may  be  observed  for  a case 
where  the  optical  path  lengths  in  the  test  section  are  known  (e.g«, 
n " constant).  Ihis  picture  is  referred  to  as  the  "no-flow"  or  reference 
interferogram  and  is  denoted  by  the  subscript  "o" . Then 


where  p ^ 

Jj  ”Po'“l'o 

is  known,  and  No  is  measured  from  the  no-flow  Interferogram. 


After  the  event  being  studied  in  the  test  section  is  recorded  by  the  "flow" 
interferogram,  may  be  used  in  Eq.  (U.3)  to  find  n-i  providing  no  changes 
have  occurred  in  the  optical  path  lengths  outside  the'test  section  (l.e., 

*’’*  ■ However,  this  assuBqpticm  is  often  subject  to  unceartainty, 

particularly  if  the  time  between  the  recording  of  the  "no-flow"  and  "flow" 
Interferograms  is  sigalficant.  A less  stringent  dependence  on  the  no-flow 
picture  can  be  invoked  by  the  use  of  a "reference  point"  in  the  flow  inter- 
ferograms as  described  later. 

Consider  the  coordinate  system  shown  in  Fig.  10  where  the  light  propaga- 
tion is  initially  in  the  z-directlon.  In  the  most  general  case,  the 
refractive  index  may  vary  over  the  entire  field,  such  that  n » n(x,y,z). 

The  light  path  must  also  be  such  that  the  travel  time  be  a minimum  as  given 
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by  Fermat’s  principle. 


•a' 


J n(x,y,z)di  = 0 


where  denotes  the  first  variation  of  the  integral. 

However y to  keep  the  analysis  from  becoming  overly  cumbersome,  the 
following  assumptions  are  made  for  a first  order  approximation: 

1.  To  satisfy  Fermat's  principle,  the  refraction  of  light  is 
small  such  that  each  beam  may  be  considered  to  travel  in 
a straight  line  path  in  the  z-direction,  perpendicular  to 
the  (x,y} -plane. 

2.  The  refractive  index  n varies  only  across  the  field  of  view 
or  in  the  (x,y} -plane,  and  does  not  vary  along  the  beam 
direction  z. 

3>  Ibe  test  section  windows  enclosing  the  test  chaniber  are 
exactly  parallel,  such  that,  for  any  (x,y) -value,  the  geo- 
matrical  distance  across  the  teat  section  in  the  z-dlrectlon 
la  the  same  value  L. 

The  first  assumption  is  satisfactory  for  many  flows  being  examined,  but 
becomes  questionable  in  regions  where  density  gradients  are  high.  The  second 
merely  emittiaslzes  that  the  intex*ferometer  is  most  convenient  in  analyzing 
two-dimensional  flows , 'and  becomes  much  more  difficult  to  use  quantitatively 
if  changes  occur  in  the  beam  direction  because  of  the  integrated  effect.  The 
third  assumption  is  usually  good,  and  may  be  pre-oontrolled  to  a known 
accuracy  in  the  equipmant. 

Mathematically,  the  above  stlpxilations  translate  into  £q.  (4.3)  in  the 
following  manner: 

1.  d/^  ■ dz 

P 

2.  n^  ■ n^(x,y),  t*  • Ta(x,y) 

3.  b-a  ■ L for  all.  (x,y) 


1 


Note  that  the  order  of  interference  N may  vary  across  the  field  of  view  or 
■become  a function  of  (x,y)  if  the  optical  path  lengths  outside  the  test 
section  and  the  refractive  index  inside  the  test  section  vary  across  (x,y). 
Equation  (li.U)  makes  it  possible  to  have  a digital  analysis. 

Some  useful  comments  concerning  the  adjustment  of  a Mach-Zebnder  inter- 
ferometer will  be  made.  In  the  "ideal”  case  with  perfect  optical  components 
and  coherent  li^t,  it  would  be  poesible  to  adjust  the  interferometer  perfectly 
such  that  the  optical  lengths  outside  the  test  section  do  not  vary  over  the 
(x,y) -plane  or  field  of  view,  such  that  t*  is  constant,  and  not  a function 
of  (x,y).  Thin,  with  constant  conditions  in  the  test  section,  there  would  be 
a constant  interference  across  the  (x,y} -field  and  hence  a uniform  illumination 
over  the  entire  viewing  screen.  In  other  words,  the  screen  would  be  totally 
bright,  dark  or  gray.  This  adjustment  leads  to  what  is  termed  the  "Infinite 
fringe"  method  (i.e.,  the  fringes  are  spaced  infinitely  far  apart). 

When  changes  take  place  in  the  test  section  during  a flow,  any  fringes, 
or  lines  of  extreme  interference,  which  are  visible  will  Just  correspond  to 
lines  of  constant  refractive  index.  This  follows  immediately  from  Eq.  (4.4). 


n(x,y)  L + tw  ■ N(x,y)  X 


Since  a constant,  the  orders  of  Interference,  N(x,y),  which  are  seen  on 
the  screen  are  directly  related  to  the  refractive  index,  n(x,y).  Along  any 
fringe,  N(x,y)  ■ constant,  and  therefore  n(x,y)  must  be  a constant.  This 
method  is  often  used  for  a qualitative  examination  of  a flow,  particularly 
in  oases  where  the  refractive  index  can  be  directly  related  to  the  density  of 
the  medium,  as  in  a perfect  gas.  Any  fringes  which  are  visible  are  Isopycnics 
or  lines  of  constant  density.  However,  the  fringes  are  usually  quite  wide 
and  accuracy  is  therefore  impaired.  Consequently,  the  infinite-fringe  method 
is  poor  for  quantitative  measurements  that  require  very  good  spatial  resolution 
of  density.  Since  only  the  extreme  lines  of  interference  (fringes)  are  easily 
detected  from  Interferograms , the  refractive  index  and  hence  density  must 
change  a sufficient  amount  through  the  flow  to  cause  changes  in  optical  path 
lengths  equivalent  to  Integer  changes  of  wavelength.  Otherwise,  perhaps 
only  one  or  two  fringes  will  be  visible  over  the  entire  flow  from  which  to 
determine  the  results.  In  other  words,  the  spatial  resolution  is  generally 
poor.  This  also  makes  it  difficult  to  determine  the  exact  order  of  inter- 
ference along  the  visible  fringe,  which  is  necessary  in  calculating  the  re- 
fractive index.  Finally,  it  is  almost  impossible  to  obtain  the  perfect  optical 
components  to  provide  vtniform  interference  over  a very  large  area.  Consequently, 
a fringe  or  two  may  appear  only  because  of  optical  Inperfections , and  this 
effect  may  be  difficult  to  account  for. 

These  problems  are  overcome  by  using  the  "virtual  fringe"  method.  In  this 
case,  one  or  more  of  the  components  of  the  interferometer  are  rotated  slightly 
to  produce  ^cajetrical  path  differences  between  the  test  section  and  ccapen- 
sating  beams  outside  the  test  section.  The  basics  of  this  technique  are  shown 
in  Fig.  11,  in  which  the  clockwise  rotation  of  splitter  83  by  an  angle  e about 
the  x-axis,  produces  different  distances  in  the  light  path  directions  to  where 
the  interfering  beams  reunite.  For  positions  (y,z)  on  the  splitter  above  the 
exis  of  rotation,  the  path  lengths  of  the  test  section  arm  are  increased  ^ile 
those  of  the  compensation  arm  are  shortened,  and  below  the  axis,  the  opposite 
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is  true.  Looking  over  the  (x^)-field  of  view,  as  y increases,  the  path 
lengths  in  the  test  section  arm  continuously  increase  linearly  while  those 
in  the  conpensation  arm  continuously  decrease,  independent  of  x.  Therefore, 
with  constant  conditions  in  the  test  section,  horizontal  lines  of  extreme 
interference  will  appear,  each  successive  fringe  occurring  at  a value  of  y 
where  the  path  length  difference  l,s  again  such  as  to  produce  the  next  minimum 
or  maximum  as  shown  on  actual  interferograms  in  Fig  .4a.  As  the  angls' 
increases,  the  path-length  changes  in  the  y-direction  become  more  seyere 
and  the  spacing  between  the  fringes  becomes  smaller.  Each  successive  fringe 
in  the  y-direction  corresponds  to  the  next  higher  order  of  interference,  since 
the  test  section  path  lengths  Increase  as  y increases . Had  the  rotation  of 
been  in  a counter-clockwise  direction,  the  changes  as  a function  of  y would 
be  reversed  and  horlzontao.  fringes  would  again,  appear,  although  the  orders 
of  interference  wotild  now  decrease  in  the  y-direction.  Although  this  very 
simplified  description  of  fringe  formation  omits  the  details  of  focussing 
(Ref.  8),  it  should  suffice  for  this  analysis. 

With  available  .controlled  rotation  of  the  optlceJ.  components  about  both 
the  X and  y-axes,  the  spacing  and  orientation  of  these  ’Virtual"  fringes  can 
be  chosen  for  the  optimum  spatial  resolution  of  Lhe  flow  being  studied.  The 
more  fringes  that  are  in  the  region  of  interest,  the  more  locations  are 
known  at  which  there  is  identifiable  Interference.  Consider  one  such  location 
(x,y)  where  there  is,  for  exasple,  after  r'»iustment  of  the  interferometer,  a 
maximum  corresponding  to  a particular  integer  order  of  interference  H‘.  If  the 
refractive  index  in  the  test  section  then  changes  in  the  region  of  (x,y),  the 
optical  path  lengths  \d.ll  change  altering  the  order  of  interference,  and  there 
will  no  longer  be  a maximum  at  this  location.  !Ihe  maximum  corresponding  to  N' 
will  now  be  at  some  new  location  where  the  new  optical  path  length  inside  the 
test  section  and  the  appropriate  path  difference  outside  again  cause  the  same 
constructive  interference.  It  is  the  measurement  from  an  Interferogram  of 
fractional  changes  in  fringe  position  which  permits  a sensitive  determination 
of  changes  in  refractive  index  and  hence  density. 

In  an  Interferogram  such  as  Fig.  6,  the  lines  of  constant  extreme  Inter- 
fennce,  or  fringes,  along  which  there  is  total  constructive  or  destructive 
interference,  are  easily  Identified.  For  example,  each  bright  line  corresponds 
to  a maximum,  where  N(x,y)  ■-5',  and  N*  is  some  integer  along  the  line. 

Althouj^  the  exact  integer  value  N ' , may  not  be  known , the  fringe  can  be 
assigned  an  arbltary  "fringe  nxiniber"  M'.  The  next  bright  line  of  interference 
will  generally  correspond  to  an  Integer  value  of  either  N'  t i or  N’  -1, 
depending  on  the  direction  of  rotation  in  the  initial  virtual  fringe  adjust- 
ment. In  other  words,  between  consecutive  bright  fringes  (or  between  consecutive 
dark  fringes),  M*  * 1.  Hence,  if  this  next  fringe  is  assigned  a fringe 

nuidjer,  M*  * 1,  such  ^at  • 1,  there  will  be  a linear  relationship 
established  between  the  assigned  fringe  numbers  M',  and  the  actual  Inteference 
orders  N’.  In  mathematical  terms,  if  N*  ■ + -6M*,  then  H’  « + M’  + k, 

^rtiere  k « constant.  ~ ” 

If  the  fringes  have  been  given  increasing  numbers  in  the  same  direction 
as  the  actual  orders  of  interference  increase,  then  the  plus  sign  prevails, 
while  if  the  numbering  were  done  in  the  opposite  direction,  the  minus  sign 
must  be  used.  If  each  fringe  has  been  assigned  an  integer  fringe  nuniber 
correctly  relative  to  the  others,  there  will  be  the  same  different  k between 
the  true  orders  of  interference  and  the  assigned  fringe  numbers  over  the 
entire  interferogram.  Thus,  the  fringe  numbers  will  represent  the  actual 
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orders  of  inteference  to  Vi thin  the  iinknown  constant  k.  It  should  he  pointed 
out  ftt  this  stage  that  the  major  problem  in  analyzing  complex  flows  from 
interferograms  is,  in  fact,  the  correct  assignment  of  fringe  numbers  over  the 
field  to  represent  the  interference  correctly.  Often,  fringes  become  difficult 
to  follow,  and  extreme  care  must  be  taken.  This  aspect  will  be  discussed 
later. 

The  fringe  number  M*  may  be  generalized  to  a continuous  linear  veuriable 
M,  as  was  done  for  the  interference  order  N,  in  Eq.  (4,3).  In  this  case, 
the  centre  of  a fringe  is  taken  as  the  location  where  M has  the  integral 
value  of  the  fringe  number.  This  corresponds  to  the  assumption  that  the 
position  of  extreme  interference  (either  total  constructive  or  total 
destructive)  is  at  the  centre  of  the  observed  brightness  or  darkness.  Then 


N(x,y)  « + M(x,y)  + k 


Substituting  into  Eq.(4.4)  yields 


n(x,y)L  + T*(x,y)  ■ + M(x,y)X  + (4.6) 


From  any  interferogram,  M(x,y)  can  be  determined  after  numbering  each 
successive  fringe  correctly.  If  the  point  (x,y)  lies  at  the  centre  of  a 
fringe,  then  M(x,y)  will  be  the  integer  fringe  number,  whereas  if  (x,y)  lies 
anywhere  between  the  centres  of  two  bounding  fringes  M(x,y)  will  be  an 
interpolated  value  between  the  two  fringe  numbers . Either  the  dark  or  bright 
fringes  may  be  numbered  consecutively,  provided  that  AM  > 1 between  a 
corresponding  type  of  fringe.  For  more  accuracy,  each  bright  fringe  could  be 
given  an  integer  value  and  each  dark  fringe  between  given  a half-integer 
value,  although  this  procedure  vould  double  the  work  of  the  an^ysis. 

The  digital  analysis  of  any  interferogram  involves  determining  M(x,y) 
at  any  desired  location  on  the  picture.  This  procedure  may  be  broken  into 
three  basic  parts  as  follows. 

Assignment  of  Fringe  Nunibers 

The  fringes  over  the  entire  interferogram  are  initially  numbered  relative 
to  each  other  as  correctly  as  possible.  When  beginning  this  process  on  any 
picture,  the  number  assigned  to  the  first  fringe  is  completely  arbitrary 
since  the  actual  order  of  inter *erence  is  unknown.  As  mentioned,  however, 
once  a particular  fringe  has  been  given  a number,  all  other  fringes  must  be 
numbered  correctly  relative  to  it  in  order  that  a constant  difference  k 
(unknown)  be  maintained  over  the  entire  Interferogram  between  the  selected 
fringe  numbers  and  the  actual  orders  of  interference.  Each  interferogram 
will  have  its  own  value  of  k which  can  be  easily  handled  in  the  subsequent 
analysis . 

Figure  12(a)  is  an  example  of  a no-flow  interferogram  taken  with 
uniform  conditions  and  hence  uniform  refractive  index  in  the  test  section. 
Consecutive  numbering  of  the  dark  fringes  was  done  in  a straightforward 
manner  as  shown,  beginning  at  the  entirely  arbitrary  value  of  21.  The 
assignment  of  fringe  numbers  for  any  no-flow  picture  is  usually  quite  routine 
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unless  optical  inperfections  are  severe. 

Figure  12('b)  is  the  subsequent  flov  interferogram  in  which  an  ionizing 
shock  wave  in  argon  is  travelling  through  the  test  section  from  right  to 
left  as  it  diffracts  around  an  expansion  corner.  (It  may  be  compared  with 
Pig.  8 where  vertical  fringes  were  used  to  reveal  different  flow  details.) 

It  is  important  to  point  out  that  the  number  assigned  to  the  first  fringe 
in  this  ](>icture  is  again  conpletely  arbitrary,  and  need  have  no  relationship 
whatsoever  to  the  no-flow  interferogram.  To  enphasize  tlis  point,  the  fringe 
nxanbers  have  been  purposely  chosen  different  from  the  no-flow  picture  even 
in  the  pre-shock  region.  It  is  advantageous,  however,  to  number  both  pictures 
in  the  same  direction.  The  necessary  criterion  is  only  that  the  relative 
numbering  between  fringes  on  a particular  interferogram  be  correct . In  this 
case,  some  care  must  be  taken  xn  identifying  and  numbering  the  fringes  near 
the  wall  and  through  the  shock  wave.  There  are  more  complex  flows  for 
which,  in  certain  regions,  the  correct  numbers  may  be  very  difficult  to 
ascertain,  and  intuition  or  experience  is  necessary ♦ However,  it  should  be 
mentioned  that  a mistake  in  numbering  of  even  one  fringe  will  usually  give 
physically  unrealistic  results.  Hence,  the  correct  integer  values  may  become 
obvious  from  the  results,  and  the  interferogram  may  be  easily  re-analyzed 
w^th  the  proper  fringe  nximbers. 

Digital  Storage  of  Fringe  Coordinates 

The  locations  of  each  fringe  of  interest  over  the  Interferogram  are 
then  recorded  using  a suitable  "digitizer”.  A digitizer  consists 
basically  of  a "tablet"  on  which  a chart  or  picture  may  be  secured,  and  a 
travelling  cursor  which  may  be  moved  at  will  over  the  tablet.  The  location 
of  the  cursor  at  any  point  on  the  tablet  is  displayed  continuously  by  a 
digital  read-out  of  the  x and  y coordinates  (usually  In  one  hundredths  of 
an  Inch)  with  respect  to  the  axes  and  origin  of  the  tablet.  These  (x,y}- 
coordinatas may  be  recorded  through  some  interface  such  as  digital  tape, 
either  continuously  at  a desired  rate  of  points  per  second,  or  point  by 
point,  by  depressing  a "log"  button  with  the  cursor  at  the  desired  location. 

The  interferogram  is  placed  on  the  digitizing  tablet,  with  its  (x,y)- 
coordinate  system  oriented  correctly  with  respect  to  the  (x,y)-axeB  of  the 
tablet . This  is  generally  accomplished  by  using  known  surfaces  or  cross- 
wire  locations  on  the  interferogram  as  a reference  (see  Figs.  8 and  12)  to 
ensure  the  correct  positioning  and  provide  the  necessary  correspondence 
between  the  digitizer  coordinate  system  and' that  used  on  the  interferogram. 
After  alignment  has  been  established,  the  centre  of  each  fringe  is  then 
traced  over  the  interferogr im  using  the  digital  cursor . When  recording 
in  the  point  mode,  the  cursor  is  moved  in  steps  along  the  fringe,  de- 
pressing the  log  button  at  each  step  to  record  the  (x,y)-coordlnate . In 
the  continuous  mode,  the  log  button  is  held  depressed  and  the  recording 
of  coordinates  is  done  automatically  at  a pre-selected  number  per  second  as 
the  cursor  is  moved  continuously  along  the  fringe.  The  two  recording  modes 
are  used  interchangeably  with  continuous  recording  in  regions  where  the 
fringe  is  smooth,  and  a more  careful  point  by  point  recording  where  rapid 
variations  occur.  In  either  mode,  the  recorded  points  must  be  sufficiently 
close  for  the  subsequent  analysis  such  that  the  fringe  location  between  any 
two  points  can  be  accurately  given  by  a Joining  straight  line.  Consequently, 
when  digitizing  has  been  coiipleted,  a sufficient  number  of  points  will  have 
been  digitized  (recorded)  to  accurately  define  the  location  of  each  fringe 


anywhere  over  the  interferogram. 

An  example  of  the  results  of  this  procedure  is  shown  in  Fig.  13.  In  this 
case,  the  dark  fringes,  digitized  from  the  interferogram  in  Fig.  12 (h)  have 
been  replotted  vising  an  IBM  l627  Calconp  Digital  Plotter  connected  with  an 
IBM  1130  conputing  system.  Each  fringe,  of  course  is  actually  described  by 
a number  of  recorded,  discrete  points  as  shown  at  the  top  cf  the  figure  for 
a particular  fringe  number.  In  the  full  plot,  the  plotter  pen  has  mere]y 
been  left  down  as  it  travelled  from  point  to  point  along  each  fringe.  As 
can  be  seen,  it  is  an  accurate  representation  of  the  interferogram,  and  the 
data  is  now  in  digital  form  suitable  for  analysis. 

Determination  of  M(x.y) 


The  value  of  M at  any  desired  location  (x,y}  is  necessary  if  Eq.  (it .6)  is 
to  be  used  to  find  the  refractive  index  at  that  point.  Since  many  points 
across  the  interferogram  have  been  recorded  by  the  digitizing  procedure  at 
which  integer  values  of  M (fringe  numbers)  are  known,  a two-dimensional 
interpolation  scheme  would  seem  appropriate,  if  not  somewhat  lengthy,  in 
finding  M at  other  positions.  However,  if  the  position  of  each  fringe  can 
be  determined  accurately,  even  between  the  discrete  points  by  linear  inter- 
polation (straight-line  segments)  as  specified  above,  the  problem  remains  only 
in  finding  M(x,y)  at  off -fringe  locations.  Restating  this,  since  lines  can 
be  defined  from  the  digitizing  procedure  along  which  M is  known,  all  that 
remains  is  to  determine  M at  positions  not  on  the  lines.  This  is  accomplished 
in  the  following  manner. 

The  interferogram  is  divided  into  a grid  of  constant  "argumant”  lines, 
as  closely  spaced  as  desired  and  generally  perpendicular  to  the  average 
direction  of  the  no-flow  fringe*.  The  most  convenient  choice  of  argument 
lines  (although  not  a necessity)  for  data  handling  are  lines  alcmg  which 
one  coordinate  remains  constant.  For  example,  with  basically  horizontal 
fringes  as  in  Fig.  13  a grid  of  lines  of  constant  "x"  is  generated  as  shown 
in  Flg.lli(a).  For  vertical  fringes,  a grid  of  lines  of  constant  y would  be 
most  suitable.  From  the  digitized  data,  the  intersection  points  of  any 
fringes  with  an  argument  line  can  be  found.  Thus,  for  each  argument  line,  a 
plot  may  be  made  of  fringe  number  versus  the  coordinate  of  intersection.  Ihe 
qtviantlty  M may  then  be  found  at  any  coordinate  along  the  argument  line  by  . 
a one-dlmenslonal  Interpolation.  As  a result,  to  find  M at  any  location 
(x,y),  the  correct  argument  line  pas'jlng  through  the  point  is  chosen,  along 
with  discrete  values  of  M are  known,  and  a one-dimensional  interpolation  is 
done  along  the  line  to  find  M at  this  point. 

The  procedure  is  demonstrated  in  Figs.  iMa)  and  (b)  wliere  the  grid 
consists  of  lines  x.  ■ constant.  For  each  argument  line  a plot  may  be  of 
the  variation  of  fringe  number  M with  y as  shown  in  Fig,  l4(b).  An  inter- 
polation of  chosen  order  may  be  done  to  determine  M at  any  specific  value 
of  y.  This  process  is  handled  easily  by  a conqputer  and  may  be  performed 
quickly  and  often. 

With  a method  available  for  determining  M(x,y)  as  desired  from  any 
interferogram,  the  use  of  Eq.  (U.6}  will  now  be  examined  in  more  detail. 

After  the  interferometer  has  been  suitably  aligned  with  the  desired  number 
and  spacing  of  virtual  fringes,  a no-flow  interferogram  is  taken  where  the 
refractive  index  in  the  test  section  has  a known  and  uniform  value. 


n(x,y)  “ Dq  “ constant 


Then 


n^L  + T^*(x,y)  = + M^(x,y)X  (*+•?) 

vhere  again » the  subscript  o denotes  the  no-flow  picture.  Although  n is 
known  and  M (x,y)  can  be  determined  from  the  interferogram,  k the  ° 
constant  difference  between  the  fringe  numbering  and  the  real^orders  of 
interference  is  not  known. 

For  the  subsequent  interferogram  taken  of  the  flow  under  study,  Eq. 
(4.6)  will  be  sinqply  restated: 


n{x,y)L  + T*(x,y)  » + M(x,y)X+  kX 

It  is  advantageous 
same  direction  as 
both  expressions, 

[n(x,y)  - n J L + 

If  no  changes  have 
two  interferograms 

is  known,  and  M(x,y)  and  M (x,y)  may  be  measured^from  the  respective  inter-*^ 
ferograms,  it  is  then  possible  to  determine  the  refractive  index  in  the  flow 
n(x,y),  to  within  the  constant,  (k-k^),  i.e,, 

[n(x,y)  - n J L ■ + [K(x,y)-M^{x,y)3X  « (k-k^)X  (4.9) 


that  the  fringe  numbering  in  the  flow  picture  be  in  the 
the  no-flow  picture , such  that  the  same  sign  holds  in 
Subtracting 

C T*(x,y)-T^*(x,y)]  » + [M(x,y)  - M^(x,y)]X  (4.8) 

+ (k-ko)  X 

taken  place  outside  the  test  section  during  the  time  when  the 
are  taken,  then  T*(x.y)  ■ Ti'*(x.y).  Moreover,  since  n 


If  the  difference  between  the  chosen  fringe  numbering  and  the  actual  orders 
of  interference  happens  to  be  the  same  in  both  pictures,  then  k ■ k , and  the 
calculations  may  be  made  directly.  If  not,  it  is  known  that  (k-k  ) Sust  be  an 
Integer . and  it  is  not  difficult  in  most  cases  to  find  its  value  from  the 
results  since  an  incorrect  choice  usually  gives  unrealistic  answers. 

However,  it  is  found  from  practical  experience  that  the  assumption  of 
no  change  in  optical  path  lengths  outside  the  test  section  between  the  no-flow 
and  flow  Interferograms  is  often  Incorrect.  This  usually  arises  from  the 
inability  to  control  the  outside  conditions  absolutely.  Small  temperature 
changes  (and  hence  density)  and  minute  ]^.8lcal  movements  of  the  optical 
components,  caused  either  thermally  or  otherwise,  are  difficult  to  avoid. 
Becavise  of  this , a stronger  method  is  often  used  which  requires  that  a 
reference  point  exist  somewhere  on  the  flow  interferogram  at  which  the  re- 
fractive index  is  known.  The  use  of  this  reference  point  can  at  least 
accomodate  uniform  changes  inT*(x,y)  between  no-flow  and  flow  pictures  and 


— 


eliiQlnates  the  ziecessity  of  choosing  the  correct  value  for 

Denoting  the  coordinates  of  this  reference  point  as  (x  ,y  ),  Eg.  (^.8) 
must  he  satisfied  at  this  particular  location  since  it  is  valid  for  any 
point  on  the  interferogram . 


[ ^ ^ *(*r*^r^  “ 1 * lLM(x^.yj.)  - M^(x^.yy) 


+ (k-k^) ^ 


Subtracting  this  expression  from  the  general  Eg.  (4.8)  for  all  point 


[n{x,y)  - n{x^,y^)]L +[T*(x,y)  - r^x.y)]  - tT*(Xj..yj.)  - ^ 


» + C[M(x,y)  - M^(x,y)]  -[M(x^,y^)-  M^(Xy,yy)  ]3X 

(4.10) 

If  the  changes  which  take  place  in  T*^(x,y)  during  the  time  the  no-flov  and 
flow  pictures  are  taken  are  uniform  for  all  (x,y),  then 

T»(x,y)  - T^(x,y)  ■ T»(Xj.^j.)  - •^o*(Xy.yy) 


(It  should  he  noted  here  that  it  is  almost  imposilble  to  account  for  changes 
which  might  take  place  outside  the  test  section  which  vary  across  the  field 
of  view.)  Eguation  (4.10)  then  reduces  t(^; 


[ n(x,y)  - M(x,y)  - M^(x,y)] 


-[M(x^.y^)  -M^(x^,y^)]I 


(4.11) 


Therefore,  using  some  reference  point  (x^,y^): 

(a)  n(x  ,y  ) must  he  known,  hut 

r r 

(h)  k>k^  need  not  he  known, 

(c)  n need  not  he  known. 

0 

The  term  M(x^,y^)  - *^Q(x^,y^)  is  a particular  value  determined  from  the 

interferograms . It  not  only  accounts  for  the  refractive  index  change  in  the 
test  section  at  the  particular  point  (x^.y^)  hut  also  for; 


(a)  Tbe  fringe  numbering  difference  'between  the  no-flov 
and  flow  pictures,  and 

(b)  any  uniform  changes  occurring  outside  the  test  section  during  the 
time  the  two  pictures  are  taken. 


This  is  easily  seen  if  the  reference  point  can  be  chosen  in  a region  tdiere 

conditions  in  the  test  section  have  not  changed  between  the  no-flow  and  flow 

plct\u*es,  as  in  the  region  ahead  of  the  shock  wave  in  Fig.  12(b).  Then,  any 

non-zero  value  of  M(x  ,y  ) - M (x  ,y  ) must  be  due  to  either  the  choice  of 

r r o r r ' * — ■ 

fringe  numbering  made  in  each  picture,  or  to  changes  outside  the  test  section. 

The  right-hand  side  of  Eq.  (^.ll)  is  termed  the  fringe  shift  S which 
coiqjhares  the  observed  change  in  interference  at  the  point  (x,y}  between  the 
no-flow  and  flew  interferograms  to  that  at  a reference  point  (x^,y^). 


S{x,y)  « + {[M(x,y)  - M^(x,y)]  - [M(x^,y^)  - M^(x^,y^)])  (U.12) 


Ad  discussed  previously,  the  choice  of  sign  depends  on  the  initial  virtual 
fringe  adjustment  of  the  Interferometer.  If  the  fringe  numbers  have  been 
arbitrarily  chosen  to  increase  in  the  same  direction  as  the  actual  orders  of 
Interference  increase,  the  plus  sign  holds,  whereas,  if  the  numbering  were 
done  in  the  opposite  Erection,  the  minus  sign  must  be  used.  If  the  Incorrect 
choice  of  sign  is  made,  all  fringe  shift  values  will  be  opposite  in  sign  from 
reality,  and  the  calculated  refractive  index  will  then  vary  from  the  reference 
value  in  the  reverse  manner.  Often,  negative  values  of  n(x,y)  result,  easily 
identifying  the  mistake.  In  almost  all  cases,  there  is  some  region  in  the 
picture  where  it  is  known  that  the  density  and  hence  refractive  index  either 
Increases  or  decreases,  such  as  thro\]gh  a shook  wave,  compression  wave, 
expansion  wave,  contact  surface,  slipstream  and  boundary  layer.  Consequently, 
the. correct  sign,  which  holds  over  the  entire  Interferogram,  can  be  easily 
confirmed  or  corrected  from  a brief  look  at  results , and  presents  no  problem 
in  the  analysis. 

In  summary,  the  usual  fringe  shift  equation  for  interferometry  may 
be  written  from  Eq.  (4.12)  as: 

P(x»y)  ■ ^ [n(x,y)  - n(x^,yy)]  (4.13) 


where 


S(x,y)  - + C[M(x,y)  - M^(x,y)]  -[M(Xy,  y^)  - M^(Xy,y^)]) 


Using  a digital  analysis  technique,  M and  M may  be  determined  from  the 
respective  no-flow  and  flow  interferograms  at  any  location  (x,y).  Including 
(x  ,y  ) thereby  allowing  a calculation  of  S.  Therefore,  knowing  L,X,  and 

enables  a final  determination  of  refractive  index,  n,  at  any 
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location  (x,y) . 
5.  CONCLUSIONS 


-rhe  experimental  techniques  associated  with  the  UTIAS  10  cm  x 18  cm 
Njrpervelocity  Shock  Tube  were  described.  The  inaccuracy  involved  in  the 
measurements  of  the  initial  pressure  P , initial  temperature  T^,  incident 
shock-wave  Mach  number  and  the  density  P at  any  point  of  a given  inter- 
ferogram  were  calculated. 

A new  computerized  method  for  evaluating  complex  flow  inter ferograms 
was  presented.  The  present  method  has  the  following  advantages  over  previous 
techniques : 

1)  It  can  easily  handle  cosplex-flow  fields. 

2)  It  can  account  for  the  changes  in  the  optical  path  length 
that  occur  outside  the  test  section  during  che  time 
interval  taken  to  record  the  "no-flow"  and  "flow"  inter- 
ferograms . 

3)  Once  the  interferograms  are  digitized,  corrected  and 
stored  on  a computer  disk,  the  total  density,  electron- 
number  density  and  degree  of  ionisation  profiles  can  be 
obtained  along  any  predetermined  line  wi-^nin.  several 
minutes . 

One  outstanding  problem,  whinh  would  reduce  the  evaluation  even  further, 
is  to  develop  a totally  automated-computerized  method  of  evaluating  Inter- 
ferograme  rather  than  the  semi -computer! zed  (hand  digitized)  method 
presented  here.  Undoubtedly,  such  a solution  would  evolve  if  more 
experimenters  made  more  use  of  Interferometry  for  investigating  complex 
flows.  However,  the  correct  assignment  of  fringe  number  over  the  field  of 
view,  to  represent  correctly  the  interference  order,  would  probably  still 
be  done  manually. 
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TABLE  1: 


IHITIAL  CONDITIONS  FOR  NGN-COMBUSTlOH  ft.. 


M 

8 

'■i.i 

Driving  Gas 

Test  G&s  > Oxygen 

1.9 

40 

2.1 

80 

°2 

2.5 

270 

°2 

2.8 

660 

3.5 

130 

hI 

U.2 

230 

He 

U.9 

670 

He 

5.7 

1610 

He 

7.2 

930 

h 

7.8 

1610 

«2 

Test  Gas  - Nitrogen 

2.0 

no 

CO 

3.7 

220 

H 2 

4.8 

690 

He 

6.2 

800 

He 

7.0 

1550 

“2 

7.8 

3320 

4 

Test  Gas  - Argon 


2.0 

60 

CO 

3.0 

70 

He 

4.4 

350 

He 

5.2 

800- 

He 

6.1 

520 

H2 

6.9 

1000 

4 

7.9 

1900 

4 

The  thickness  of  the  loylar  diaphragm  is  calculate 


TABLE  2:  GLADSTOHE-DALE  CONSTAIITa  FOR  AROOM  AND  KRYPTON 


constant 


Z, 


Argon 

(Kj«  0.65  K^) 

0.157i» 

0.1629 

3.314 

0.8718 

Krypton 

0.1127 

0.1188 

1.599 

0.4244 

(Kj-  0.75  Kj^) 


1)  X 

2)  Xg  “ 3»+71.5  S 


^1*  ^2*  \ ^2  om'^/g. 


TABLE  3 ; GLADSTONE-DALE  CONSTANTS  FOR  OXYGEN  AND  NITROGEN 


constant 


Oxygen 


\ \ \ \ h h 

11  2 2 ^ ^ 


0.1907  0.185  0.1988  0.192  0.0057  O.OO68 

+ 1!<  +1% 


Nitrogen  0.2376  0.328  0.2k6  0.331 

+ 0.03<  + 1.7IJ  + 0.03Jf  + 2.5 


1)  - 6943  % 

2)  X2-3it7l.5X 

1^,  Kg^,  ^1  ^2  in  cm^/g. 
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FIG.  1 GENERAL  TOP  AND  SIDE  VIEW  OF  THE  UTIAS  10  CM  X 18  CM  HYPERVELOCITY  SFK)CK  TUBE. 


Travelling  Time  fram  0-t(/ii9ec) 
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PIG.  3 INSTRUMENTATION  SET-UP  FOR  RECORDING  THE  NON- STATIONARY 
PHENOMENA  IN  THE  TEST  SECTION. 


FIG.  4 DISTANCE-TIME  DIAGRAM  OF  THE  INCIDENT  SHOCK-WAVE  AS  IT  MOVES 
ALONG  THE  SHOCK  TUBE.  TIMES  ARE  MEASURED  FROM  D (see  Fig.  3). 

THE  SIX  POIn'TS  CORRESPOND  TO  STATIONS  D,F,G,H,I  AND  J OF  FIG.  3. 
THE  SOLID  LINE  REPRESENTS  THE  BEST  FIT  ACROSS  THESE  MEASilREMENTS. 
THE  VELOCITY  OF  THE  INCIDENT  SHOCK-WAVE  AT  ANY  POINT  IS 


Measured  Initial  Mixture  Pressure  (H2  ,02  , He  ) = 423  psi  (Heise  Gauge) 
Calculated  Diaphragm  Bursting  Pressure  = 2452  psi 


] I* — 5 msec 


a)  DRIVER  PRESSURE  at  CLOSED  END 


Measured  Shock  Mach  Number  * 15.0  at  Station  M 
Measured  Initial  Pressure  Pe~6.l2torr  (0.12  psi) 

Predicted  Equilibrium 
Preeture  ■ 38.0  pel 
(1965  torr)  ^ 

Predicted  Frozen 
Pressure  ■ 33.3  psi 
( 1722  torr ) 


Po 

— I [* — !0  yiisec 

b)  SHOCK- WAVE  PRESSURE  at  STATION  M 
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FIG.  5 TYPICAL  PRESSURE  HISTORIES  AT  THE  TEST  SECTION  (a)  AND  IHE  BACK  END  OF  THE  DRIVER 
SECTION  (b). 


FIG.  6 SCHEMATIC  OF  MACH-ZEHNDER  INTERFEROMETER 
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FIG.  7 SCHEMATIC  DIAGRAM  OF  THE  LIGOT  SOURCE  AND  CAMERA  OPTICS  FOR  OVERALL  MAGNIFICATIONS  OF 
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FIG.  9 BASIC  COMPONENTS  OF  A MACH-ZEHNDER  INTERFEROMETER. 


FIG.  10  COORDINATE  SYSTEM  USED  FOR  INTERFEROMETRIC  ANALYSIS 


QMHM.  POaTION  IV 
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FIG.  11  'VIRTUAL'  FRINGE  CREATED  IN  TOE  PLANE  OF  SPLITTER  S 


FIG.  14  (a)  SUPERPOSITION  OF  ARGUMENT  LINES  x*  - CONSTANT  ON  DIGITIZED 

INTERFEROGRAMS.  (NOTE  IN  AN  ACTUAL  ANALYSIS  THE  GRID  LINES 
MOULD  BE  MUCH  MORE  CLOSELY  SPACED) . 


(b)  PLOTS  OF  FRINGE  NUMBER  M VERSUS  FRINGE  POSITION  y ALONG  ARGU- 
MENT LINES  xj. 
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U.  SUPPLEMENTARY  NOTES 
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20.  abstract  (Continue  on  roverao  tide  If  ncc.sJAf^  mnd  Identify  by  block  number) 

The  UTIAS  10  cm  x 18  cm  Hypervelocity  Shock  Tube  has  been  used  in  recent  years 
to  study  ionizing  shock  structures,  flat-plate  and  side-wall  boundary  layers  and 
nonstationary  oblique  shock-wave  diffractions  over  compression  corners.  These 
phenomena  were  recorded  using  a 23-cm  dia  field  of  view  Mach-Zehnder  interfero- 
meter equipped  with  a giant-pulse  dual-frequency  ruby  laser.  In  order  to  extract 
the  maximum  amount  of  data  in  these  complex  flows  a digital  evaluation  method 
was  employed.  For  this  technique  a new  approach  to  the  theory  of  interference 
was  developed.  In  this  approach  the  spatial  coordinates  (x,y)  of  the  various 
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lines  cf  interference  (fringes)  on  the  interferograms  are  put  into  digits'*  form, 
th'^reby  making  a computer  analysis  possible.  exp^<iriment«l  te^luiiquafar^ 

in^  trumentation  associated  with  the  various  measurerients  involved  Tir  rls6drctt* 
using  the  UTIAS  10  cm  x 18  cm  Hypervelocity  Shock  Tube  are  described.  Finally 
the  maximum  possible  absolute  and  relative  errors  associated  with  these  measure- 
ments are  calculated  for  four  different  gases  (argon,  krypton,  oxygen  and 
nitrogen)  which  are  usually  used  as  test  gases. 
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